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ANALYTICAL STUDY TO DEFINE A HELICOPTER 


STABILITY DERIVATIVE EXTRACTION METHOD* 


By John A. Molusis 
Sikorsky Aircraft Division 
United Aircraft Corporation 


SUMMARY 


A method is developed for extracting six degree-of- freedom stability and 
control derivatives from helicopter flight data. Different combinations of 
filtering and a priori derivative estimate are investigated and used with a 
Bayesian approach for derivative identification. The combination of filtering 
and a priori estimate found to yield the most accurate time response match to 
flight test data is determined and applied to CH-53A and CH-54B flight data. 
The method found to be most accurate consists of (l) filtering flight test 
data with a digital filter followed by an extended Kalman filter (2) identify- 
ing an a priori derivative estimate with a least square estimator and 
(3) obtaining derivatives with the Bayesian derivative extraction method. 

The results demonstrate (l) the importance of the a priori estimate, 

(2) the least square estimator used with the optimal data filter provides 
a time history match nearly as well as the Bayesian estimator, ( 3 ) identified 
derivatives only approximate conventional quasi-static derivatives and, ( 4 ) the 
identified derivatives yield the correct characteristic roots. 


INTRODUCTION 


The ability to identify aerodynamic stability and control derivatives 
from aircraft flight data has long been recognized as a desirable goal. 
Improvement in AFCS design, quantification of stability and control character- 
istics, and improvement of analytic prediction methods are some of the 
potential benefits to be realized. Identification of potential stability 
problems in prototype aircraft and flying qualities assessment can also 
benefit from derivative identification. Because of the complexity and coupled 
behavior of the helicopter, there is strong motivation for identifying the 
aerodynamic derivatives that describe flight behavior. For example, present 
analytic prediction methods show discrepancies in predicting the high speed 
unstable dutch roll roots associated with the CH-53A helicopter. Additionally, 
high gains of rate feedback used in the AFCS cause rotor tip path plane 
oscillations that are not predicted accurately by present analytic methods. 
Stability derivative identification can provide the potential tool for 
isolating the causes of these observed discrepancies. 

* The contract research effort which has led to the results in this report was 
financially supported by USAAMRDL (Langley Directorate). 
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Many applications of derivative identification from fixed wing and VTOL 
aircraft have appeared in the literature. 1 * 2 >3 These studies have generally 
been characterized by three degree-of- freedom models, negligible input driving 
noise (process noise), and sufficiently long data records with properly excited 
modes. The helicopter and many VTOL aircraft do not fall into this relatively 
simple characterization. In particular, the typical helicopter requires a 
minimum of six degrees of freedom and contains significant process noise. 
Finally, long data records often are unavailable for proper identification. 

Many applications of derivative identification have concentrated on curve 
fitting an assumed mathematical model to test data by minimizing a quadratic 
function of the fit error. Gradient techniques, least squares methods, and 
Newton Raphson methods have been used to accomplish the required minimization. 
The resulting solution generally provides a good fit to the data; however, the 
identified value of the derivative may depart considerably from the real value 
of the derivative. This problem can be attributed to improper modelling of 
the problem, insufficient data, and multiple solutions due to the algorithm 
used. 


Current statistical methods of identification provide the means to treat 
all the problems normally associated with derivative extraction from test data. 
The Bayesian approach to estimation treats the derivative identification 
problem as a stochastic problem. The criterion to be satisfied is to determine 
the most probable estimates of the stability derivatives and state variables, 
given noisy measurements. This formulation has advantages in the identification 
problem. First, we are attempting to determine the most probable derivatives 
and not simply solve a curve fit problem. Secondly, the approximations made 
in obtaining a computationally practical algorithm can be justified, or errors 
resulting from these approximations minimized. Finally, derivative identifi- 
cation from flight test data is indeed a statistical problem, and problems of 
insufficient data and proper mode excitation can be interpreted conveniently. 

An excellent review of methods for system identification can be found in 
Reference 4. In particular, the Bayesian method is developed and a recursive 
algorithm provided for the solution. This method is seen to be identical with 
the extended Kalman filter when the measurement equations are linear. This 
method is also called a maximum a posteriori estimator and is related to the 
classical maximum likelihood estimator by Bayes* rule. Reference 4 discusses 
this relationship in more detail. Reference 5 provides applications of a 
modified Kalman filter algorithm for VTOL aircraft derivative identification. 

The results demonstrate the superiority of the extended Kalman filter method 
over standard curve fit methods of Newton Raphson and conjugate gradient. 

The helicopter derivative identification problem using Kalman filtering 
methods is described in detail in Reference 6. This paper explored the 
particular aspects of helicopter derivative extraction that require special 
consideration and proposed methods for their solution. The paper demonstrates 
(l) the need for more than one maneuver when using short data records, (2) the 
existence of significant levels of process noise, and (3) the non-unique 
character of the six degree -of -freedom derivative model. In addition, the 
Bayesian maximum likelihood method is extended to process any number of 
simultaneous maneuvers , and a method for obtaining an a_ priori derivative guess 
is presented. 
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This report presents applications of the methods of Reference 6 to actual 
helicopter flight test data- An overall procedure, including filtering and 
a_ priori derivative guess , is developed to yield a helicopter derivative 
extraction method. The method is applied to CH-53A and CH-5^B flight test 
data and demonstrates the requirements necessary for successful identification. 


SYMBOLS 


A 
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A* 


B° 


B* 


Matrix of partial derivatives obtained by partially 
differentiating f(X, Xp) with respect to X. 

Matrix of partial derivatives obtained by partially 
differentiating f(X, Xp) with respect to Xp. 

Matrix of stability derivatives used to initialize 
the Bayesian max. likelihood estimator. 

Matrix of stability derivatives obtained from the 
Bayesian max. likelihood estimator. 

Matrix of control derivatives used to initialize 
the Bayesian max. likelihood estimator. 

Matrix of control derivatives obtained from the 
Bayesian max. likelihood estimator. 


A ls> B ls Lateral and longitudinal cyclic control input, 
respectively, degrees at the head. 


f(X> Xp) 


Functional relationship among the state variables and 
derivatives in the stability derivative model. 


* 1 > *1 


Lower cutoff frequency used for determining noise statistics 
from the data, Hz. 


f 2 , F 2 


Upper cutoff frequency used for determining noise statistics from 
the data, Hz. 


f“c Actual cutoff frequency used on the data, Hz. 

G Process noise gain matrix. 

ll> l2> Distance to the center of gravity of the airflow measurement 
I3 probe. Ft. 

L Normalized rolling moment defined in equation (l). 

Lp Stability derivative 3L/3P defined in equation (2). 
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M 


MB ls 

Mp 

Mq 

n 

N 

Nq 

P 

*11 

P lp 

Ppp 

P°, P* 


P 

q 

Q 

r 

R 


t 

t o» tf 

U ’ V x 

v, V 

y 

v 

w(t) 

w, v. 


Normalized pitching moment defined in equation (l). 

Stability derivative 8M/3B]_ S defined in equation (2). 

Stability derivative 3M/3p defined in equation (2). 

Stability derivative 3M/3q defined in equation (2). 

Zero mean, white gaussian measurement noise. 

Normalized yawing moment defined in equation Cl). 

Stability derivative 3N/3q defined in equation (2). 

Covariance matrix for the Kalman filter. 

State variable covariance matrix for the MLE. 

Cross covariance matrix for the state and parameter vector of 
the MLE . 

Parameter covariance matrix for the MLE. 

Initial and final parameter covariance matricies for the MLE. 
Aircraft roll rate, radians/sec. 

Aircraft pitch rate, radians/sec. 

Process noise intensity matrix. 

Aircraft yaw rate, radians/sec. 

Measurement noise intensity matrix. 

Time, Secs. 

Initial and final time respectively, secs. 

Aircraft longitudinal velocity. Ft. /Sec. 

Aircraft lateral velocity, Ft. /Sec. 

True Airspeed, Ft. /Sec. 

Zero mean, white gaussian process noise 

Aircraft vertical velocity. Ft. /Sec. (positive up) 

Longitudinal normalized force defined in equation (l). 


k 


X System state vector for the Kalman Filter and MLE. 

Xp Vector made up of stability derivatives. 

Y Lateral normalized force defined in equation (l). 

Z Vertical normalized force defined in equation (l). 

Z Measurement variable. 

Zw Stability derivative 3Z/9W defined in equation (2). 

a > 8 Aircraft angle of attack and sideslip respectively, degrees. 

<5 Perturbation from trim condition. 

<f>» ©> V Aircraft roll, pitch and yaw attitudes respectively, radians. 

Positive direction is right roll, up pitch and right yaw. 

0 coll. Collective control input, degrees at the head. 

0 tail Tail rotor collective pitch, degrees. 

Subscripts: 


( )i 
( ) T 
( r 
(*) 


Refers to the ith maneuver number. 

Denotes matrix transpose. 

Denotes matrix inversion. 

Denotes time derivative, d( ) 

dt 


Miscellaneous : 

All stability and control derivatives are presented in either 
normalized or dimensional form. The units are defined in 
Table 8 and Table 9* Dimensional derivatives are obtained by 
multiplying normalized derivatives by the appropriate mass or 
inertias . 

Roll Rate I, Pitch Rate I and Yaw Rate I used in the time 
history figures refer to the inertial aircraft rates <j>, 0, 
respectively. 

The data supplied by USAAMRDL (Langley Directorate) was chosen 
at a 1+5 knot condition and consisted of h maneuvers. Subsequently 
a one maneuver case was selected which is trimmed at 32 knots. All 
titles for this data, however, indicate a 1+5 knot trim. 
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PROBLEM FORMULATION AND METHOD FOR SOLUTION 


The primary difficulty in obtaining accurate stability derivatives from 
flight test data is due to the noise contaminations present. Additionally, 
helicopter six-degree-of-freedom derivative models are not uniquely defined 
due to the lumped effects of the rotor. This section will discuss these 
particular problems of helicopter derivative identification, present the deri- 
vative model which is used throughout this report and discuss in detail the 
Bayesian approach and procedure used. 


Review of System Identification Requirements 

Helicopter stability derivative extraction requirements differ from 
fixed wing or VTOL aircraft primarily in the following ways: 

1. The six degree-of- freedom helicopter derivative model consists of 
the lumped effects of body and rotor. 

2. The plant driving noise and measurement noise are larger for 
helicopters. 

3. The number of degrees of freedom to describe helicopter rigid 
body response is generally greater. 

h. The helicopter has unstable phugoid roots and on occasion unstable 
dutch roll roots . 

Reference 6 and Reference 7 describe some of the above differences in 
detail. A summary of the conclusions and methods for treating these problems 
will be discussed. 

Derivative definition . - The helicopter when interpreted as a six 
degree-of- freedom derivative model does not possess a unique set of derivatives 
during actual flight conditions . This is because the rotor effects are assumed 
to be lumped with the body derivative and since the rotor is always in motion 
the lumped derivative reflects a combined body and 'average 1 rotor contribution. 
The conventional method for obtaining helicopter 6 degree-of- freedom derivatives 
is based on the assumption that a perturbation in a body attitude, rate or 
velocity causes the rotor tip path plane to reach a steady state trim. The 
lumped effect of rotor and body causes a force or moment which divided by the 
perturbation yields the quasi-static derivative. 

It should be clear that under normal flight conditions (pilot excitation 
and wind gust always taking place) the helicopter does not actually fly like 
the quasi-static model. This is because the rotor is never in a steady state 
condition. The six degree-of- freedom model which ’best” represents the flight 
behavior is what is identified from flight test. The identified derivatives 
can be thought of as an ’’average” or ’’nominal” derivative. The quasi— static 
and derivative identified from flight data are further compared below. 
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Quasi-Static _ Body + Rotor Contribution after tip 

Derivative Contribution path plane reaches steady state 


System Identified = Body + Average Rotor 

Deri vat i ve Cont r ibut i on Co nt r ibut i on 

Reference 6 demonstrates that the two derivative models yield nearly the 
same characteristic roots and both models yield good time history reproduction 
of the non-linear helicopter simulation from which they were obtained. 

Correlation of quasi-static derivatives obtained from a non linear com- 
puter simulation against six degree-of-freedom derivatives identified from 
flight data can not be expected to yield identical results. However , meaningful 
characteristic root comparisons can be made as well as time response comparisons. 
Additionally, correlation of derivatives identified from flight test and com- 
puter models can be made as long as the same derivative identification method 
is used. 

The fact that six degree-of-freedom derivative models are not uniquely 
defined is primarily due to the choice of a six degree-of-freedom model rather 
than a larger more accurate representation. It is possible to approximate 
quasi-static derivatives from flight data but requires the identification of 
a nine degree-of-freedom model which includes body and tip path plane degree 
of freedom. The identified nine degree-of-freedom model can then be perturbed 
just as a non-linear computer simulation and quasi-static derivatives obtained. 

The studies conducted in this report are concerned only with six degree- 
of-freedom identification. The resulting derivative model can be used to obtain 
characteristic roots and time history generation. 

Noise contamination . - There are several types of noise present in heli- 
copter test data which are of particular concern to the derivative identification 
problem. These are as follows: 

1. Measurement noise resulting from the inaccuracy of the physical 
transducer used, 

2. Measurement Bias resulting from steady drift or calibration errors, 

3. Process noise resulting from all effects which physically drive the 
equations of motion (included are higher degrees of freedom which 
are not modeled, non linearities which are not modeled, wind gust 
disturbances and servo actuator errors), 

k. Trim errors resulting from the inability of the pilot to perfectly 
trim the vehicle in flight. 

Measurement noise can be accounted for in the problem formulation, however 
the larger the uncertainty in the measurement the longer the data record required 
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in the identification. This is why filtering the data first improves the 
derivative estimates for a fixed data length. 

Measurement bias can he formulated as part of the identification problem 
however this increases the size of the state vector and increases computation 
time. Alternately the bias can be identified first by using an extended Kalman 
filter. This approach has the advantage of providing good estimates of the 
bias as well as excellent estimates of the system state variables. This latter 
approach is the one used in this study. 

Process noise can also be formulated as part of the identification 
formulation Process noise has two effects on the identification. First, the 
larger the process noise the more data and input excitation required. Second, 
process noise tends to make the linearizations used in the derivative identifi- 
cation algorithm less accurate. Additionally process noise is represented in 
the formulation by white noise with gaussian statistics. This requirement may 
not always be met , however by using long data records or many short data 
segments this requirement should be generally satisfied. 

Trim errors can be included in the formulation however at the expense 
of increased computations • For this reason it is important for the vehicle to 
trim accurately or to fly about a trim condition sufficiently long so that 
the trim can be easily identified by visually inspecting the data. This is 
particularly important when using several independent segments of data 
(multi -maneuver case) in the identification. 

Number of degrees of freedom . - As indicated earlier the helicopter 
requires a minimum of six degrees of freedom. Present algorithms for deriva- 
tive identification are capable of treating this size problem without 
excessive computation time. The nine degree-of- freedom problem is dimension- 
ally too large for most methods including the Bayesian technique used in this 
report. However a slight modification to the Bayesian approach can provide 
a practical algorithm for treating almost any number of degrees of freedom. 

This can be done by identifying one row of derivatives at a time while 
accounting for the errors due to the remaining derivatives by a bias correction 
term. This approach will not be discussed any further here. This method 
however can be used to identify very large derivative arrays. 

Unstable characteristic roots . - The helicopter possesses unstable 
phugoid roots and on occassion has unstable dutch roll roots at high speeds. 
This does not cause problems in the identification, however when correlating 
the identified derivative model against test data, the error between the time 
history match is governed by an unstable error equation. The roots of this 
error equation are the same as the roots of the identified derivative model. 
Thus , the smallest amount of process noise or initial condition error will 
cause the linear model response to diverge from the test data even if the 
correct derivatives are identified. For this reason application of identifi- 
cation techniques will generally provide a much better time response match for 
fixed wing vehicles which generally have stable characteristic roots. 
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The Derivative Model 


The six degree-of-freedom derivative model used throughout this study 
is given below: 
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Y = 

Z = 

L = analogous to equation (2) 

M = 

N = 


It should be noted that the inertial terms which are usually found on 
the left of the equal sign in equation (l) have been lumped into the derivative 
terms. The derivatives are thus total force derivatives which are the sum of 
the aerodynamic and inertial derivatives. 

The derivative model is conveniently represented in state variable form 
by defining the state vector and parameter vector below: 
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The derivative model of equation (l) then reduces to the state vector 
equation 

X = f(X, Xp) (k) 

Equation (4) is the derivative model used in all the Bayesian maximum 
likelihood identification studies. 


Method For Derivative Identification 
(The Bayesian Max. Likelihood Estimator) 

The derivation of the Bayesian maximum likelihood estimator for simulta- 
neous maneuvers is given in Reference 6. The initials MLE are used throughout 
this report to represent the Bayesian maximum likelihood estimator and not the 
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classical maximum likelihood estimator. The actual formulation and equations 
used will he presented here. 

The Bayesian Max. Likelihood Estimator determines the most probable 
estimates of the state variables and stability derivatives given the measured 
data. The equations relating the noise contaminations, stability derivative 
model and measurements are presented below s 


System equations . 

x.(t) = f(x., Xp) + G. w. (t), Xi(to) 

h = 0 Xp(t 0 ) 

where, i - represents the maneuver number. 

Measurement equations . 

Z± = Xi + n. 


(5) 


( 6 ) 


The best estimate of the state variables and stability derivatives is 
obtained from the equations below. 


The Bayesian Max. Likelihood Estimator . 
Filter Equations : 

X = f(X, Xp) + PppR -1 ( Z - X ) 

4 = P pl R _1 (Z - x ) 


( 7 ) 


Covariance Equations : 

p n = A n p n + ( A 11 P ^) T + A ip p pi + ( A ip p pi) T - PnR - 1 P-n + Q 


■ n" r n 


P pl * P plAll T ♦ P ppA lp T - 

p = _p R-lpT 


( 8 ) 


pp ■ pi* r pi 

The individual terms in equations (j) & (8) are defined in Appendix A. 
Initialization of The MLE and Noise Statistics . 

Equations (7) and (8) require the following quantities before a solution 
can be obtained; 

(1) Measurement Noise Covariance - R 

(2) Process Noise Covariance - Q 
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(3) ' Initial State Estimate - X (t Q ) 

(M Initial State Covariance - P 1:L (t 0 ) 


(5) Initial Derivative Estimate - X (t 0 ) 



Physically the measurement noise covariance represents the amplitude of 
the power spectrum of the assumed white noise contamination on the state 
variables . Figure 1 shows a spectrum of white noise and the rigid body signals 
which the derivative model is to reproduce. 


The calculation of the measurement noise covariance R can be automated 
by filtering the raw data at frequencies f-^ and f 2 and computing the 
corresponding mean square value of the filtered data. The measurement noise 
covariance is then 


R = 


( q 2 2 - g l 2 ) 

* (f2 - fi) 


(9) 


where ; 


°2 = the mean square value of the raw data filtered at f 

2 

= the mean square value of the raw data filtered at 


The process noise covariance can be computed in an analogous manner by 
operating on the accelerations . Since the accelerations contain both process 
noise and sensor noise, the sensor noise should be subtracted out. A reasonable 
estimate is obtained without subtracting the sensor noise and this is the 
approach used in this study to obtain process noise covariance. 

The initial conditions on the system equations were obtained from the 
first data point of the measured data. The initial state covariance is then 

p ll(t 0 ) = Rirf c (10) 

where; f c = the filtering cutoff used on the data. 

R = computed from equation ( 9 ) for each state. 

The initial derivative guess and variance in the guess was obtained using 
two methods. First, a least square estimator was employed which operates on 
data filtered by an extended Kalman filter. The least square estimator provides 
an estimate of the derivatives and variance. The variance required modification 
in each case. The least square estimator is discussed in Appendix B. The 
second method employes engineering judgement to select the derivative guess 
and variance. 
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The Data Filters 


Three filtering methods are investigated in this study as follows: 

1. First order low pass with cutoff freq. of 10 cps (8 cps for the 
CH-54B data). 

2. Graham Digital Filter with cutoff frequency (FC) of 3 cps and 
termination frequency (FT) of 4 cps. 

3. A Kalman filter (state variable estimator) which minimizes the 
overall uncertainty in the data and removes bias error. 

The first order low pass filter with cutoff frequency at 10 cps was 
chosen to represent the case in which little or no filtering is done to the 
data. Since all instrumentation contains inherent filtering and an aliasing 
filter is needed before discretizing the data, the 10 cps filtered data re- 
presents essentially no filtering with respect to rigid body modes (typically 
less than U cps). 

The Graham digital filter represents a high quality zero phase shift 
filter which has extremely high roll off rate. The cutoff frequency of 3 cps 
and termination frequency of k cps essentially retains only rigid body 
frequencies and low frequency contaminations. 

The Kalman filter is the highest quality filter of the three in that it 
accounts for random discrepancies of related measurements and also can remove 
bias error. The Kalman filter is really a state variable estimator. The 
predicted state estimates do not in general filter out high frequency data. 
Many of the high frequency effects are a part of the true motion of the 
vehicle due to higher modes of motion and it is only the engineering modelling 
requirement to obtain rigid body data that the high frequency data is labeled 
as undesirable. Since rigid body motion is what is desired , all data used by 
the Kalman filter is first filtered with the Graham zero phase shift digital 
filter. Thus only low frequency data is predicted by the Kalman estimator. 

The equations describing the three filters are presented in Appendix C. 

Determination of the Noise Statistics for the Kalman Filter . - The 
Kalman filter requires (l) measurement noise statistics, (2) process noise 
statistics and (3) an initial state estimate and variance. Table 1 
summarizes the method used to obtain the measurement noise statistics while 
Table 2 summarizes the method used to obtain process noise statistics. 
Statistics computed from the data without engineer interaction and statistics 
inputed by the engineer are presented. Table 3 presents the equations used to 
compute the initial state and variance from the data. 



Procedure To Determine A Helicopter Derivative Extraction Method 


The procedure to determine the most accurate helicopter derivative 
extraction method is shown in the flow diagram of Figure 2. Three different 
filtering techniques are employed; a first order low pass, a Graham digital and 
a Kalman filter. Additionally two different methods for obtaining the a priori 
derivative guess and variance are used; a least square method and by using 
engineering judgement. Combinations of the above filtering methods and 
a priori guess are used in conjunction with the Bayesian Max. Likelihood 
Estimator. A discussion of the criteria used to determine the most accurate 
method is presented in a subsequent section. 


TEST DATA USED AND HELICOPTER DESCRIPTION 


Flight test data from two different helicopters Eire used in this study; 
a CH-53A at both 100 knots and 150 knots trim condition and flight data from 
a’ CH-54B at 45 knots trim which was supplied by AAMRDL, Langley Research Center; 
The most accurate method is determined from CH- 53 A data at 100 knot trim. 


The data used in this study is summarized below. 


Flight Test 
Data 


No. of Independent 
Maneuvers used 


Comments 


CH-53A 

(100 kts , 352"C.G.) 


6 (6 secs, length each) 

(4 maneuvers in Identification, 
2 maneuvers used to determine 
the most accurate method) 


CH-53A 

(150 kts, 352"C.G.) 4 (5 seconds length each) 


CH-54B 4 (6 seconds length each) 

(45 kts, 346" C.G.) 


CH-54B 1 (16 seconds length) 

(45 kts, 346" C.G.) 


This case was used to 
determine the most 
accurate method. 


using most accurate 
method 

using most accurate 
method 

using most accurate 
method 


CH-53A Helicopter 

The flight test data from the CH-53A was obtained from a handling 
qualities flight test program conducted for the Navy by Sikorsky Aircraft. 
Extensive dynamic stability data are recorded on magnetic tape and selected 
maneuvers from the data are used for derivative identification. 


20 


Description of Test Vehicle . - The test aircraft is a CH-53A helicopter 
modified into a CH-53D configuration by the installation of YT64-GE-12 engines 
and an up-rated transmission. Figure 3 shows a photo of a CH-53A in flight. 
Table 4 gives the physical characteristics of the test CH-53A helicopter. 

Test Data at 100 Knots Using the Three Filters . - Six maneuvers are 
selected with at least one of the primary control inputs excited for each 
maneuver. All four control inputs are excited in the four maneuvers used in 
the identification. The inputs used are typically 1 degree - \ second 
duration pulse inputs and thus a data record length of 6 seconds ' duration is 
used to prevent large deviations from linearity. 

Figure 4 through Figure 9 show the six maneuvers of CH-53A data at 
100 knot trim condition filtered at 10 HZ. with the first order low pass filter. 
These data are not corrected for c.g. location and represent the outputs of the 
physical transducer. These same data filtered with the Graham digital filter 
at 3.0 HZ can be found in Volume II (Figure 1 through Figure 6). 

Figure 10 through Figure 15 show the same six maneuvers of CH-53A data 
(100 kts) converted to body axis. The Kalman and digital filtered data are 
shown superimposed. Examination of these figures reveal the ability of the 
Kalman filter to remove bias error and make use of the coupled information in 
several channels to predict a more accurate estimate. 

Test Data at 150 Knots Using The Three Filters. - Figure 1 6 through 
Figure 19 show the four maneuvers of CH-53A flight data at a 150 knot trim 
condition filtered ..at__lQ__ Hz jwith the first order . low. pass filter. These 
data are not corrected f or _c.g. location and represent the outputs of the 
physical transducer. These same data filtered with the Graham digital filter 
at 3.0 Hz can be found in Volume II (Figure 7 through ^Figure 10.). 


Figure 20 through Figure 23 show the same four maneuvers of CH-53A 
data (150 knots) converted to the body axis. The Kalman and digital filtered 
data are shown superimposed. The Kalman filter is shown with bias error 
removed and the uncertainty due to random errors minimized. 


CH-54B Helicopter 

The flight test data from the CH-54B was supplied by Langley Directorate , 
USAAMRDL , and stored on magnetic tape. Selected maneuvers from this tape are 
used to identify derivatives. Two cases are under investigation; identification 
of derivatives from four 6-second simultaneous maneuvers and identification 
from one l6-second maneuver. 

Description of test vehicle . - The test aircraft is a CH-54B Army crane. 
Figure 24 shows a photo of a CH-54B in flight. Table 5 gives the physical 
characteristics of the test CH-54B helicopter. 
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Test data at k nots using the three filters (U maneuver case) . — Pour 
maneuvers are selected which contain excitation from al I four input controls . 

The four maneuvers chosen contained in some cases dramatic differences in trim 

condition. These maneuvers were selected primarily because all control inputs 
were excited and good excitation was not available in other segments of data 
that was supplied by Langley Directorate, USAAMRDL . The large differences in 
the trim condition led to a poor identification as was later accessed by resim- 
ulation of the identified model and comparing against the test data. Since this 
case contained improper trim all test data and results are presented in Volume 
II of this report. 

Volume II (Figure 11 through Figure l4) show the CH-54B flight test data 
(four maneuver case) filtered at 8 HZ. with the first order low pass filter. 

These data are not corrected for c.g. location and represent the outputs of 
the physical transducer. These same data filtered with the Graham digital 
filter at 3.0 HZ. is shown in Volume II (Figure 15 through Figure 18). 

Volume II (Figure 19 through Figure 22) shows the same four maneuvers of 
CH-54B data (45 kts.) converted to the hody axis. The Kalman and digital 
filtered data are shown superimposed. 

Test data at 45 knots using the three filters (one l6— second maneuver 
casej_. - A single l6-second data record from the CH-54B flight data was selected. 
All control inputs are excited and a well defined trim was apparent from the 
data. Figure 25 shows these data filtered at 8 HZ . with the first order low 
pass filter , uncorrected for c.g. location and represent the outputs of the 
physical transducer. These same data filtered with the Graham digital filter 
at 3.0 HZ. is shown in Volume II (Figure 23). 

Figure 26 shows the l6-second maneuver of the CH-54B data (45 kts.) 
converted to the body axis. The Kalman and digital filtered data are shown 
superimposed# The Kalman filtered data is shown with bias error removed and 
the large random error appearing in the linear velocities is also shown 
removed . 

DETERMINATION OF THE MOST ACCURATE HELICOPTER DERIVATIVE EXTRACTION METHOD 


Ten methods consisting of different data filtering and a priori derivative 
guess were investigated in this study. The purpose of this investigation is to 
examine the three data filters previously discussed, access the least square 
method and determine the best method for obtaining an ji priori derivative 
estimate for initializing the maximum likelihood method. Table 6 lists the 
ten methods used in determining the best overall approach. The arbitrary 
derivative guess noted in Table 6 was obtained as follows . The control deriva- 
tives were obtained by examination of the test data, measuring the accelerations 
where the pulse inputs reached their peak and dividing the accelerations by 
the magnitude of the control input pulse. The stability derivatives Zq and Y r 
are known approximately and have a magnitude of 169 ft / s ec . (approximately the 
trim condition.) All other stability derivatives were set to zero. The 
variance on the derivatives were guessed using engineering judgment. Table 7 
summarizes the derivative value and variance for the "arbitrary" guess . 
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Criteria for Selection of the Most Accurate Method 


The most accurate derivative identification method is determined by using 
four maneuvers of CH-53A flight test data at 100 knot trim condition in the 
identification process. The identified derivative model is then simulated and 
compared against two independent maneuvers not used in the identification. 

The root mean square (RMS) error between the time histories of the identified 
derivative model and test data are then computed. RMS error for both test 
data used in the identification and not used in the identification is computed. 
The RMS error for all state variables and acceleration time responses are 
summed and the method yielding the lowest resultant RMS error for the two 
maneuvers not used in the identification is selected as the optimum method. 

The resultant RMS error for the four maneuvers used in the identification is 
used as a secondary check for selecting the optimum method. 

Results Using the Ten Methods 

The ten methods of Table 6 were used to obtain stability derivatives 
from the four maneuvers of CH-53A flight test data previously presented. 

Table 8 presents the values of the identified stability derivatives while 
Table 9 presents the identified control derivatives. The time history compari- 
sons of the ten methods investigated are presented in Volume II with the ex- 
ception of the method found to be most accurate. The ten methods are pre- 
sented in the following figures: 


Method 

Volume 

Figure 

1 & 7 

ii 

2U-29 

2 & 8 

ii 

30-35 

3 & 9 

ii 

36-Ul 

4 

ii 

42-47 

5 

ii 

1+8-53 

6 

ii 

54-59 

3 & 10 

i 

27-32 


The Most Accurate Method 

The most accurate method was chosen by computing the RMS error between 
the test data (Kalman filtered) and the output of the identified derivative 
model. This was done for data used in the identification and data not used in 
the identification. The results sire shown in Table 10. Table 10 shows that 
Method 10 yields the smallest RMS error for both data used in the identifica- 
tion and data not used in the identifications. 

Table 10 clearly indicates the improvement made with the Kalman Filter 
when using the least square estimator (Method 1, 2 and 3). Also the Kalman 
Filtered data gives superior results when using the MLE with the arbitrary 
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initial derivative guess (Method 4, 5 and 6). Contrary to the above, the 
Kalman Filter did not yield the best result when using the MLE with the least 
square initial derivative guess. The best result according to the RMS criteria 
used was obtained with the digital filtered data into the MLE. This conflicting 
fact is attributed to two factors. First, the RMS error criterion used does 
not include the effects of wind gust disturbances and thus can not be expected 
to yield consistant results when these disturbances are present. Second, the 
small RMS error difference between Method 9 and Method 10 requires more data 
to provide statistically consistent results. 

It is felt at this stage of the study that the method outlined above 
(Method 10 ) is the best approach using the 100 knot data and detailed results 
of this method are presented below. However, application of both Method 9 and 
Method 10 to the 150 knot data shows the opposite conclusion; that is Method 
9 provides the best result. This will be discussed under the section 
Application to CH-53A at 150 knots . 

The RMS error criteria can not be expected to yield consistant results in 
the presence of wind gust and for unstable models. This method of evaluation 
is particularly unsatisfactory for comparing small differences in RMS error. 
Since both Method 9 and Method 10 yield very close RMS error values the most 
accurate method could not be clearly identified. At this stage of the study 
both Method 9 and Method 10 are acceptable even though Method 10 shows a 
slightly improved RMS error. Method 10 is selected for the 100 knot CH-53A 
data . 


The results obtained thus lead to the following procedure for the most 
accurate helicopter derivative identification method. 

1. Process the raw data with the Graham digital filter to provide data 
with zero phase shift and frequency content only pertinent to rigid 
body response. 

2. Process the digital filtered data with the Kalman filter to remove 
bias error and provide consistent estimate of the state variables. 

3. Use these data to obtain the initial derivative guess by the least 
square estimator (LSE). Modify the derivatives' variance from the 
LSE by multiplying the variance by a constant factor (multiplica- 
tion by 100 yields good results ) . 

1+. Use the LSE derivatives and modified variance in conjunction with 
the Digital or Kalman filtered data to obtain optimal derivatives 
from the Bayesian Maximum Likelihood Estimator (MLE). 

One conclusion which is apparent from this study is that the least square 
estimator (LSE) when used with the Kalman data yields an RMS error nearly as 
good as Method 9 or Method 10. In fact, the LSE with Kalman data is superior 
to the MLE with the first order filter. In view of the fact that the LSE re- 
quires considerably less computation for solution than the MLE, the LSE may 
prove to be valuable in studies where computer time is important and particularly 
when large derivative arrays are required. 
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IDENTIFICATION OF DERIVATIVES FROM TEST DATA USING THE MOST ACCURATE METHOD 


This section discusses in detail the results obtained using the method 
found most accurate as determined in the last section. Applications to CH-53A 
data at both 100 knots and 150 knots trim and to the CH-5UB data which was 
supplied by the Army (AAMRDL) are presented. 


Application To CH-53A at 100 Knots 


The preceeding section discussed the selection of the optimum method and 
it was found that Method 10 gave the best results. This method is used to ob- 
tain the results of this section. 


Interpretation of helicopter 6 degree— of- freedom derivatives . - Before 
proceeding with the evaluation of the identified derivatives it is important 
to understand precisely how to interpret helicopter derivatives identified from 
flight test data. 

The helicopter 6 degree-of-freedom derivative represents the lumped sum 
of the rotor contribution and body contribution. The helicopter in actual 
flight conditions consists of low frequency modes of motion ( body D.O.F.) 
and high frequency modes of motion (rotor D.O.F. ). If we choose to represent 
the helicopter by a 6 D.O.F. model obtained by identification from input/output 
data, then the identified coefficients represent the lumped contribution of 
body plus average rotor. This is a different interpretation than the convention- 
al quasi-static derivative which represents the lumped contribution of body plus 
the rotor in a steady state trim condition. This difference arises because the 
helicopter/rotor has characteristic motion consisting of more than six degrees of 
freedom, and the 6 D.O.F. model is a low order approximation to actual flight 
behavior. The identified 6 D.O.F. model is a better representation than the 
quasi-static model. The helicopter rotor in flight, being continually excited, 
is not characterized correctly by static data and thus the body motion is 
more accurately approximated by the derivative made up of lumped body and 
average rotor contribution. 

The helicopter 6 D.O.F. derivatives can best be interpreted by referring 
to Figure 33. This figure was obtained from a nonlinear computer helicopter 
simulation. Figure 33 shows one typical stability derivative obtained from 
the nonlinear helicopter computer simulation. This derivative, M , is ob- 
tained by peturbing roll rate (Ap) holding all body variable constant and 
allowing the rotor to reach a new steady state trim. The resulting pitch 
acceleration divided by Ap is the conventional quasi-static derivative 
(i.e. Mp = Aq/Ap). A 6 DOF model obtained by using a system identification 
method is also shown superimposed on Figure 33. It is seen that the 6 DOF 
identified value represents the body plus average rotor contribution. This 
behavior occured with nearly all derivatives investigated and thus we conclude 
that the 6 DOF identified value is a very meaningful quantity although it 
must be compared to the quasi-static value in the manner outlined above. 
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Also superimposed on Figure 33 is an identified 9 DOF model derivative value. 
Perturbing this model results in the dashed line and reduces to a 6 DOF quasi- 
static derivative value. This figure shows that while the identified 6 D.O.F. 
derivative lends itself to physical interpretation, it requires at least a 
9 DOF model to identify the complete derivative behavior of a helicopter. 

The suggested procedure for comparing derivatives identified from flight 
data with quasi-static derivatives is to multiply the flight identified deriva- 
tive by a correction factor. The correction factor can be determined from a 
nonlinear computer simulation set up to duplicate the flight conditions of the 
vehicle under study. This correction is an approximate method to account for 
the absence of the rotor degrees of freedom. The correction is as shown 
below; 


Computer Model 
Quasi -static 
Derivative 
Computer Model * 

Identified 
Derivative 

*The computer model Identified Derivative can be approximated by 
computing the average value over one rotor revolution; of time from 
each derivative time plot (a sample plot is illustrated in Figure 33). 


Quasi-static 

Derivative 


Flight 

= Identified 
Derivative 


The identified derivatives . - Table 11 shows the values of the identified 
derivatives from CH-53A flight data (100 knots) using Method 10. These deri- 
vatives were presented earlier in Table 8 and Table 9 and are repeated here in 
a more convenient form. Inspection of the derivative array of Table 11 shows 
that Mq was identified with a positive sign. This fact is attributed to the 
a_ priori derivative value used. Derivative identification from noisy data in 
general yields nonunique solutions. A reasonably close a priori value must be 
provided to guarantee convergence to the correct solution. 


The derivative values identified in Table 11 represent the body contri- 
bution plus average rotor contribution as discussed previously. This is 
further exemplified in Figure 3^ which shows two of the identified derivatives 
(Lp and Zw) compared against theoretically predicted values. Again the complete 
derivative variation with time due to a perturbation is shown. The theoretical 
quasi-static value and 6 DOF identified value is compared. This figure 
demonstrates quite good correlation. Figure 35 shows Mp compared in an 
analogous manner and is shown identified with the opposite sign. 


Derivative convergence . - One way to access if sufficient data has 
been used in the identification is to inspect the convergence plots of all 
stability and control derivatives. Figure 36 shows the derivative Lp obtained 
using the Bayesian estimator as data is processed from one data point to the 
next . Four maneuvers of 6 seconds ' duration each were processed simultaneously 
in the identification. Derivative convergence for Lp, Mp^ and Mq are shown 
in Figure 36, Figure 37 and Figure 38 respectively. Two different a priori 
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derivative values and variance were used in each case; a least square guess 
and arbitrarily selected values . These curves demonstrate that the amount of 
data used was sufficient to allow convergence and that the better the initial 
guess and smaller variance, the faster the convergence. Reference 7 discusses 
further the implications of derivative convergence. 

Characteristic roots. - Although the derivatives obtained from a 6 D.O.F. 
identification must be interpreted differently than the quasi-static derivatives, 
the characteristic roots should be nearly the same for both identified and 
quasi— static derivatives. This assumes the effect of nonlinearities does not 
cause the time history responses to differ significantly from linearity. Figure 
39 shows the characteristic roots obtained from the least square and Bayesian 
estimated derivatives. This figure demonstrates that the least square method 
when used with Kalman filtered data yields reasonably good derivative est ima tes ■ 
Figure 1*0 shows the characteristic roots obtained from the Bayesian identified 
derivatives compared against those obtained from Sikorsky Aircraft's nonlinear 
helicopter computer simulation. Correlation is excellent except for the spiral 
mode and the lateral-to-longitudinal coupling roots. The presence of the high 
frequency rotor D.O.F. always effect the lateral-to-longitudinal roots, thus we 
conclude the spiral mode is the only mode which does not show good correlation. 


Time response comparison . - The most accurate method for identifying 
derivatives from the CH-53A at the 100 knot trim condition was discussed in 
detail previously. This consisted of a least square guess (Method 3) and the 
Bayesian estimator using digital filtered data (Method 10). The time history 
comparisons were presented in Figure 27 through Figure 33. Inspection of 
these figures reveals that generally the time response comparisons of the 
identified derivative model shows excellent fit to the primary channels 
associated with each input (i.e. B^ input shows good fit to longitudinal 
D.O.F.). The cross response (i.e. ^ input producing lateral response) 
generally shows a fair to poor match.) This behavior is due to the fact that 
only a 6 D.O.F. model is being compared. Generally it requires a higher D.O.F. 
linear model to yield good correlation to the cross response. 


Application To CH-53A At 150 Knots 

The method found most accurate (Method 10) was applied to the CH-53A 
data at 150 knot trim condition. The initial derivative guess was obtained 
from Method 3 (Least square estimator using Kalman data). The computed RMS 
error showed considerably lower RMS error for the LSE than for the MLE. (Method 
10) The LSE (Method 3) and the MLE (Method 10 ) was found to have RMSvalues of 
25.6 and 31.7 respectively. These results are in conflict with the study to 
determine the most accurate method in the last section. It was previously 
found that the RMS error for Method 10 was lower than the LSE (Method 3) which 
was used as the a priori estimate. This was true for both data used in the 
identification and for data not used in the identification. 

Method 9 was then applied to the 15 0 knot data and the RMS error was found 
to be 22.7 which is considerably less than both Method 3 (LSE) and Method 10. 

It was concluded that Method 9 provides the most accurate derivative model for 
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the 150 knot data. Figure 1+6 through Figure 1+9 shows the time response compari- 
son for Method 9. Volume II (Figures 60 through 63 ) shows the time response 
comparison of the identified derivative model (Method 2 and Method 10 ) against 
the flight test data. Table 12 compares the RMS error in the test data for the 
three identified derivative models (Method 3, Method 10, and Method 9 ). 

The identified derivatives . - Method 9 was applied to the four 5-second 
maneuvers (discussed in the section on test data used) of CH-53A data at 150 
knots trim condition. Table 13 shows the identified derivative values for the 
LSE while Table lk shows the derivative values using the MLE. The derivatives 
x u> M q> ani ^ Y v were identified with positive values and are expected to all 
haVe negative signs. This problem was discussed in the section on identification 
from the CH-53A at 100 knots from flight data and the positive sign is a result 
of not vising a close enough a priori derivative guess. 


Derivative convergence . - Figure 1+1 shows the identified derivative L 
obtained from the MLE (Method 9)- The least square guess is shown and the^ 
convergence indicates more data could have been used in the identification. 

The lateral cyclic inputs (roll inputs) are shown applied in the middle and 
also at the end of the data record. It is clear that the application of the 
input causes rapid reduction of t he uncertainty. It is emphasized that 
derivative convergence plots for all derivatives should be inspected to ensure 
enough data has been used in the identification. 

Characteristic roots. - Figure 1+2 shows the characteristic roots obtained 
from the identified derivatives for both the LSE (Method 3) and the MLE 
(Method 9)* Again the characteristic roots from the LSE are quite close to 
the MLE roots . Figure 1+3 shows the characteristic roots identified using the 
MLE from flight data compared against roots obtained from Sikorsky Aircrafts 
nonlinear helicopter computer simulation of the CH-53A at 150 knots . The 
correlation is poor and further examination is necessary. Figure 1+3 also 
shows the frequency of the dutch roll roots obtained by visual inspection of 
the flight data. This can only be done with the pitch AFCS on since the 
vehicle's motion diverges too rapidly to access the dutch roll because the 
phugoid roots are unstable. It was found that the presence of the pitch AFCS 
increases the dutch roll frequency slightly. Excellent agreement of the dutch 
roll roots from the identified derivatives compared with roots obtained by 
visual inspection of the test data (AFCS on) is obtained. 

An investigation was made to determine which derivatives are the cause 
of the unstable dutch roll roots. Figure 1+1+ shows the characteristic roots of 
the identified derivatives when the longitudinal and lateral degrees of free- 
dom are completely decoupled. It can be observed from Figure 1+1+ that the 
dutch roll roots are now stable indicating coupling as the cause of the 
instability. Also notice that the phugoid roots are now aperiodic. Thus it 
is concluded that coupling is the cause of the unstable dutch roll roots and 
oscillatory phugoid roots which is not predicted using Sikorsky Aircraft's 
helicopter computer simulation. 
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Figure b5 shows the characteristic roots of the complete coupled deri- 
vative array and the characteristic roots obtained when the derivative N q is 
set to zero. The dutch roll roots migrate toward the left half plane. Many 
other derivatives and derivative pairs were set to zero and the derivative Nq 
was found to have the most influence on the dutch roll instability. 

Time response comparison. - The four maneuvers used in the identification 
were discussed in a previous section. Figure k6 through Figure k9 show the 
time history comparison of the identified derivative models using Method 3 
and Method 9 against the CH-53A, 150 knot flight data. The time history match 
is considered a good fit. Discrepancies are apparent, however they are due 
mainly to (l) the larger number of degrees of freedom not modeled and (2) the 
highly unstable roots which cause the resimulated derivative model to diverge 
from the data even with perfectly identified derivatives. This second fact is 
discussed in detail in Appendix B. 


Application To CH-5^B At i+5 Knots 


Using four maneuvers of 6 seconds 1 duration each . - Method 9 is considered 
the most accurate method for obtaining helicopter derivatives. Method 9 is 
selected for identifying derivatives from the CH-5^B data because the RMS error 
obtained for the 150 knot CH-53A data was approximately k0% lower than that of 
Method 10, whereas for the 100 knot case Method 10 was only slightly better 
than Method 9- This method was applied to four maneuvers of CH-54B data 
supplied by AAMRDL Langley Research Center* The data was found to yield 
considerably different trim conditions and the identified derivative model 
obtained was considered a poor estimate of the characteristic motion due to 
the inconsistant trims selected. These maneuvers were selected primarily 
because all control inputs were excited and good excitation could not be 
found in any other segments of data. Volume II (Figure 6b through Figure 67) 
shows the time response comparison for the four maneuver case. Table 15 and 
Table 1 6 shows the identified derivatives using the LSE (Method 3) and the 
MLE (Method 9)* Since the four maneuver case represents an unacceptable run 
due to the poor trims, a second case consisting of one long data record was 
made. 

Using one maneuver of 1 6 seconds duration (Arbitrary selection of 
derivative variance ) . - Method 9 was applied to one long data record of CH-5^+B 
data which was discussed previously. Method 9 consists of (l) Kalman filtering 
the data, (2) obtaining an initial derivative estimate with the LSE, (3) modi- 
fying the derivative variance obtained from the LSE and {h) inputing this into 
the MLE to obtain final derivative estimates. Step (3) above was accomplished 
by arbitrarily multiplying all derivative variances by a constant factor (a 
value of 100 was found to yield the best results ) . Since in general the best 
variance to use is dependent on the particular problem under study , a different 
approach was used for selecting this variance. This was obtained by observing 
the derivatives and variance from the LSE and using engineering judgement to 
appropriately modify the derivative variance. (This method is called Method 11). 
Table 17 shows the identified derivatives using the LSE. Table l8 shows the 
1 sigma uncertainty obtained from the least square estimator after multipli- 
cation by 10 and the value selected using engineering judgement. 


29 



Method 11 was used with the derivative variance selected by engineering 
judgement (Table 18) to obtain MLE derivatives. These derivatives are presented 
in Table 19 and the time history response comparison is presented in Volume II 
(Figure 68). This figure requires further elaboration. Because the character- 
istic roots of the identified model are unstable the identified derivative model 
when compared to the test data diverges in proportion to the unstable roots . 
Figure 68 (Volume II ) thus is reinitialized after each 4 seconds because of 
this divergence. This divergence is discussed in detail in Appendix D. The 
time history comparison shows that the MLE derivative model yields a larger 
fit error than the LSE derivative model. This is a direct result of using too 
large an initial derivative variance. This problem is corrected by reducing 
the magnitude of the derivative variance. Selecting too large a derivative 
variance causes the MLE to effectively ignor the _a priori derivative estimate 
and weight heavily the data being processed. Thus a longer data record is 
needed to cause proper derivative convergence . 

Using one maneuver of 16 seconds duration (Optimum selection of derivative 
variance). - Examination of the initial derivative variance to con- 
elusion that too large a value was used. A second run was made with derivative 
variances divided by four (This is called Method 12). The resulting MLE deri- 
vatives are shown in Table 20. The time history comparison of the LSE model, 

MLE model and test data is shown in Figure 50. It is clear from this figure 
that the MLE model is now superior to the LSE model. 

Derivative convergence plots were obtained for all 60 stability and 
control derivatives using Method 12. Figure 51 presents the complete set of 
derivative convergence plots. It can be seen from Figure 51 that most of the 
derivatives have converged to within an acceptable engineering variance. 

Examination of the derivatives yields the conclusion that the collective 
and tail rotor control derivatives are different than what analytic prediction 
would show. This difference is a direct result of linearly dependent control 
inputs. Since the tail rotor input is coupled to collective stick through the 
mixing box, collective inputs cause proportional tail rotor inputs which are 
linearly dependent. The tail rotor was not excited independently with 
sufficient magnitude to remove this linear dependence. Linear dependent in- 
puts do not in general cause errors in the other derivatives or effect the 
time rssponse. 

The characteristic roots from the identified derivative model (Method 12) 
were obtained. The roots are presented in Figure 52 and demonstrate a highly 
unstable phugoid pair and an unstable spiral mode. 
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RECOMMENDED PROCEDURE FOR OBTAINING 
STABILITY DERIVATIVES FROM FLIGHT DATA 


The results obtained to date in helicopter derivative identification in- 
dicate that j. light test data should meet specified requirements for accurate 
identification. This conclusion is based upon the experience gained from 
applications to computer simulations, CH-53A flight data at 100 and 150 knot 
trim conditions and CH— 5^-B flight data. A set of guidelines and recommended 
procedures is given below which should be followed in order to obtain data 
which is suitable for 6 degree-of-freedom derivative identification. 

Preflight Considerations 

Instrumentation . - The instrumentation which should be included in the 
flight test vehicle is listed below. 

Attitude Gyros; Roll, Pitch, Yaw. 

s 

Rate Gyros; Roll, Pitch, Yaw 

Angular Accelerometers; Roll, Pitch, Yaw 

Airflow Measurements ; Airspeed, Angle of Attack, Sideslip 

Linear Accelerometers; Longitudinal, Lateral, Vertical 

Control Inputs; Longitudinal Cyclic, Lateral Cyclic, Main Rotor Collec- 
tive, Tail Rotor Collective. 

The control input measurements should he made as close to the actual 
aerodynamic surface as is possible. Primary control servo or secondary servos 
are most desirable with stick inputs considered acceptable but less desirable. 
Compound or unconventional vehicles will require different and additional 
control input measurements. 

The identification of derivatives from the hover regime would benefit 
with the inclusion of instrumentation designed to measure rotor downwash vector 
or other suitable airspeed information and is recommended where possible. 

Redundant instrumentation will always improve the identification accuracy 
and is recommended whenever the added weight and expense is justified. Redund^- 
ant linear accelerometers and/or airflow instrumentation would be the most 
beneficial in terms of improved accuracy. 

Sampling rate .- Sufficient sampling rate should be employed to avoid 
sampling errors (aliasing) and to allow for identification of rigid body short 
period dynamics. Fifty to one hundred samples per second is recommended. This 
sample rate is required to keep the numerical computation in the algorithm 
accurate. 

Filtering of data before sampling .- Standard hardware filters (Butter- 
worth , 3rcT~order~Towp^^ used on the data before sampling. 
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A high, cutoff frequency ia desirable, since, rt eliminates excessive phase shift 
near rigid body modes. Typically 10 cps to 25 cps filters are recommended. 
Additionally , a high cutoff frequency allovs automatic determination of noise 
statistics from the data as discussed in this report. 

Data storage .- Magnetic tape is the most convenient format for data 
storage. This is particularly true for processing large amounts of flight data. 

Wind gusts .- The presence of wind gusts degrades the accuracy in the 
identified derivatives. Generally the larger the magnitude of wind gust pre- 
sent, the longer the data records required to average out their effects to 
achieve the accuracy attainable without gusts. For this reason it is recom- 
mended that the flights take place during days of negligible gust or when this 
is not possible sufficiently long data records should be used to properly 
identify the derivatives . 


Flight Testing Procedure 

There are four key requirements that the pilot of the test aircraft must 
meet in order to obtain proper data for identification; (l) obtain good initial 
trims, (2) exercise the input controls to excite all modes; (3) obtain data 
records long enough to yield sufficient information for accurate identifica- 
tion and (1»J keep the vehicles motion resonably close to the initial trim to 
stay within a linear, small perturbation region. Each of the above four re- 
quirements is necessary to obtain accurate derivatives. These requirements 
will be discussed in detail with the aid of Figure 53 which gives instructions 
as to typical pilot inputs which will yield data suitable for identification. 

Aircraft trims .- Since the stability derivative model represents the 
characteristic motion of the test aircraft about a trim condition, it is im- 
portant to properly trim the aircraft before the start of a maneuver. At 
least 5 seconds duration of stick input excitation should be made to clearly 
identify a trim condition (the longer the better). Figure 53 shows this 
trimming time segment at the start of the data record for a single maneuver 
identification. When using discrete segments of data (i.e. the multimaneuver 
identification) the same initial trims for each maneuver must be obtained by 
the pilot. This includes stick positions as well as V, 8, a ,<j> and 0. 

Pilot inputs .- After the initial trim period the first pilot input c an 
be conducted. Figure 53 shows typical good inputs for identification corres- 
ponding to time t = 0 for the longitudinal cyclic. While this input is ex- 
ercised all other controls should be held at their trim. The longitudinal 
cyclic is a doublet which is positive for second then switched to negative 
for \ second and then returned to trim for 1 second. The next several. Seconds 
(2 to 5 seconds) are used to bring the aircraft back to the original trim. 

All inputs are exercised in accomplishing this. This trim should be approxi- 
mately the same as the starting trim (i.e. all stick positions and V, o, 8, 

<P >0). The doublet input allows (!) good excitation of the modes while the "return 
to trim" portion ensures the small perturbation assumption is not violated, (2) 
provides additional excitation of the modes of the system and (3) retrims the air- 
craft for the next input (Lateral Cyclic). The above procedure (i.e. doublet 
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then retrim) is then repeated for lateral cyclic, rudder and collective inputs 
respectively as shown in Figure 53. Thus all controls have been exercised both 
individually and collectively for approximately 20 seconds. The aircraft is 
now at its original trim and the pilot should now fly the aircraft about this 
trim, exercising all control inputs as shown in Figure 53. Switching inputs 
(i.e. sharp edge pulse type inputs) generally provide the best information for 
identification and these should be used where possible when retrimming and 
flying about a trim. 

Figure 53 is intended to be only a guide for generating flight test data. 
The time scale shown may not be realistic for all vehicles and thus only re- 
presents an example. The switching type inputs also represent an ideal input 
which is only an approximation. 

Length of Maneuver . - The maneuver of Figure 53 is considered, heuristic- 
ally, the most desirable and the pilot should use this figure as a guide to 
generating data for derivative identification. Figure 53 shows about 30 
seconds of data used. This should yield sufficient data for identification as 
shown by the studies in this report. Alternately discrete segments could be 
used in the identification but this requires excellent initial trims for each 
individual maneuver. 

Factors Effecting Accuracy 

Three important influences which adversely effect the accuracy of the 
identified derivatives are l) nonlinearities, 2) poor quality data and, 3) im- 
proper a priori derivative estimate. The nature of these influences are as 
follows : 

1. Nonlinearities are not modeled as part of the identification and thus 
they tend to bias some of the identified derivatives from their 
expected value. It is thus important to keep the aircraft near the 
trim condition to minimize the nonlinear effect. 

2. Poor quality data caused by either poor instrumentation or signifi- 
cant wind gust can effect the accuracy of the identified derivatives . 
Every effort should be made to use as good instrumentation as is 
considered practical and flights should be conducted during days of 
calm air. When these conditions are not met, the use of long data 
records can offset to some degree these effects. 

3. An a priori derivative estimate and variance is required to initialize 
the Bayesian maximum likelihood derivative estimation method used in 
this report. The a priori derivative variance must be chosen large 
enough so that false convergence does not occur but not so large that 
long data records are required for convergence to an accurate identi- 
fication. A good a priori derivative estimate is considered important 
when only short data records are available or when other influences 
such as nonlinearities and poor quality data as discussed above 
deteriorate the identification. A good a priori derivative guess will 
tend to offset the bias introduced in the identified derivatives due 
to nonlinearities and poor quality data. 
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An accurate a priori guess is very important for obtaining accurate 
derivatives. Since derivative identification from noisy data in 
general has a nonunique solution, a close guess must be provided to 
guarantee convergence to the correct solution. This is often true 
with the pitch derivatives Mq and Mw. Often Mq would take on a 
positive value and My a negative valve which is just the opposite 
from what analytical prediction would show. This uniqueness problem 
is eliminated by providing a accurate initial derivative guess and 
uncertainty in the guess . 

Interpretation of Identified Derivatives 

Helicopter 6 degree-of-freedom (DOF) derivatives identified from flight 
data are found to differ from conventional quasi-static derivatives . This 
difference is because the helicopter has additional rotor degrees of freedom 
which cause the rotor/vehicle to £ly differently than a 6 DOF quasi-static 
model would predict . The conclusions obtained regarding interpretation of 
derivatives identified from flight data sire summarized as follows: 

1. Helicopter 6 DOF derivatives as defined in this report represent the 
lumped contribution of the body plus rotor. The quasi-static 
derivatives assume that the rotor instantaneously follows the rigid 
body motion and, thus, the rotor contribution neglects rotor dynamic 
effects. System identified derivatives are found to consist of the 
time averaged rotor contribution over one rotor revolution and thus 
includes an average rotor contribution. 

2. Correlation of derivatives identified from flight data vs. analytic 
computer generated derivatives must take into account the consider- 
ations discussed in l) above. If a nonlinear computer model is used, 
correlation with flight identified derivatives can be made directly 
from the nonlinear computer model, since derivatives which consist 

of the body contribution and the time averaged rotor contribution can 
be obtained from the computer model. 

3. Six DOF identified derivatives yield, characteristic roots which can 
be directly compared against roots from quasi-static derivatives 
(i.e. the phugoid, dutch roll, spiral and longitudinal drag character- 
istic roots). The coupled longitudinal-to-lateral roots are generally 
affected by rotor dynamics and thus cannot be directly compared. 


Procedure To Identify Derivatives 

The derivative identification procedure found to be most accurate in this 
study consists of a) filtering the data with a Kalman filter formulated to re- 
move bias error and minimize random error while not identifying derivatives, 

b) using least square identified derivatives for the a priori estimate and 

c) obtaining final derivative estimates using the Bayesian maximum likelihood 
estimator (MLE) - This procedure is shown in block diagram form in Figure 5^. 



The overall procedure which should he followed when attempting 6 DOF derivative 
identification from helicopters is outlined below. 

1. When flight testing helicopters for derivative identification, the 
preflight considerations and flight testing procedure discussed in 
the preceeding sections should be followed. Clearly defined initial 
trims, continuously switched inputs (primarily pulse type), long data 
records and small deviations from the selected trim condition are 
considered important to successful identification. 

2. The data should be filtered with the Kalman filter to remove bias error 
and reproduce all state variables with minimum uncertainty. 

3. An a priori derivative guess should be obtained using the least square 
estimator. This guess should be modified as required. (Modification 
will be required if nonlinearity is present in the data or if poor 
quality data is obtained due to poor instrumentation or due to the 
presence of significant wind gusts.) The variance as computed by the 
LSE is always computed too small and should either be multiplied by 

a constant factor or selected using engineering judgement. 

k. Derivatives can be obtained from the Bayesian MLE using the Kalman 
filtered data and the a priori derivative estimate and variance 
determined in step 3. Derivative convergence from the MLE should be 
plotted and inspected to ensure sufficient data lengths have been 
used to yield the required accuracy in the derivatives. If the re- 
quired accuracy has not been achieved the final derivative estimates 
and variance may be used to reinitialize the MLE and a second pass 
over the same data (or new data) can be made. 

5. After derivative convergence has occurred the identified derivative 
model should be re simulated against both the data used in the identi- 
fication and data that was not used in the identification (if avail- 
able) to test the validity of the identified model. This time history 
comparison is usually valid for only b to 6 seconds length due to the 
unstable characteristic roots associated with most helicopters. The 
time history match for the six degree-of- freedom model compared 
against the test data may show error in some cases. This error is 
often attributed to the rotor degrees of freedom and unstable charac- 
teristic roots and is generally not attributed to the identified 
derivatives . 


CONCLUSIONS 


Analysis and applications of the Bayesian approach to identification of 
stability derivatives from helicopter flight data have led to the following 
conclusions : 
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i. Six, aegree-of -freedom stability and control derivative identifica- 
tion from flight test data is a feasible method for obtaining and 
evaluating the stability characteristics of helicopters. 

'd. Six degree-of-freedom derivatives identified from flight data were 
found to represent the combined body plus time averaged rotor con- 
tribution over one rotor revolution. This is an approx im ate repre- 
sentation and differs from the conventional quasi-static derivative 
interpretation, which assumes that the rotor instantaneously follows 
the rigid body. The system identified and quasi-static derivatives 
both yield nearly the same low frequency characteristic roots. 
Differences appear only in the high frequency coupled longitudinal 
to lateral roots. 

3. The procedure that yields the most accurate derivatives is as 

follows: a) filter the data with a zero phase shift digital filter 

to remove high frequency content followed by a Kalman filter to re- 
move bias error and obtain consistency of data, b) obtain an a priori 
derivative estimate and variance with a least squares estimator, 
c) modify this derivative estimate and variance to reflect a priori 
knowledge more accurately, and d) obtain final derivative estimates 
using the Bayesian maximum likelihood method. 

k. The least square estimator, when used with the Kalman filtered data, 
was shown to yield a time history match nearly as accurate as the 
Bayesian method described in 3) above. Since the least square estima- 
tor is computationally much more efficient than the Bayesium method, 
this approach deserves consideration for those problems where 
computer time is important. 

5. Factors that adversely influence the accuracy of the identified 

derivatives are : nonlinearities , poor initial trims , poor quality 

data, incorrect a priori derivative estimate and variance in the 
estimate, wind gust disturbances, insufficient data lengths, and 
improper pilot control input excitation. Each of these factors can 

be overcome or minimized to achieve accurate derivative identification. 

6. An accurate a priori derivative estimate must be used to ensure con- 
vergence to the correct derivative value. Failure to provide an 
accurate a^ priori estimate will result in identified derivatives 
which do not represent the correct physical derivative value, however, 
the characteristic roots in general should provide the correct modes 
of the system. 

7. Characteristic roots and time history response provide a useful test 
of the accuracy of the identified derivatives. Errors in the time 
response match are caused by higher degrees of freedom, wind gust 
disturbance, initial condition errors, control input trim errors, and 
errors in the identified derivatives. An identified derivative model 
with unstable characteristic roots will diverge in the time response 
match even with perfectly identified derivatives. 
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8. Results obtained to date in helicopter derivative identification from 
flight data indicate that this approach -will play an increasing role 
in understanding the coupling effects and stability characteristics 
of present helicopters. 


The list of specific conclusions presented below are a result of stability 

derivative identification from the CH-53A and CH-5hB helicopters (AF.CS off) . 

CH-53A 100 kts 

1. The CH-53A at 100 knots trim and aft C.G. shows stable dutch roll, 
unstable phugoid and a slightly unstable spiral divergent mode as 
determined from the identified derivatives. Inspection of independ- 
ent test data reveals the above characteristic roots to be correctly 
identified. 

2. Correlation of characteristic roots identified from flight data vs. 
roots obtained from Sikorsky Aircraft's General Helicopter simulation 
program shows excellent agreement except for the spiral mode . 

3. The identification of derivatives from CH-53A data confirms the 
lumped interpretation (body plus average rotor), shows generally good 
agreement against theoretically predicted values for those derivatives 
compared and yields what are believed to correct characteristic roots. 

CH-53A 150 kts 

1. The CH-53A at 150 knots trim and aft C.G. shows unstable dutch roll, 
unstable phugoid and stable spiral mode as determined from the identi- ' 
fied derivatives. Visual inspection of independent test data confirms 
the unstable dutch roll roots. 

2. Discrepancies are found to exist between characteristic roots identi- i 
fied from flight data vs. roots obtained from analytic prediction. i 
Analytic prediction shows aperiodic phugoid and stable dutch roll 
roots while identification from flight data reveals oscillatory 
phugoid and unstable dutch roll. 

_ 

CH-54B k5 kts ! 


1. Identification of one long data record was found to yield better re- 
sults than four short segments. This was primarily due to the 
different initial trims for the four segments. 

2. The identified derivatives were found to yield unstable phugoid and 
unstable spiral mode . 

3. Time history comparison of the derivative model against test data 
shows good results and clearly shows the unstable error divergence 
due to the unstable phugoid roots . 
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RECOMMENDATIONS 


The procedure developed in this report for six degree-of-freedom 
stability derivative identification has had only limited actual application. 
The potential achievable from helicopter derivative identification remains to 
be realized. The following areas of application can provide for cost re- 
ductions , better aircraft designs , and further understanding of the stability 
characteristics of present day helicopters. 

1. Use identified derivatives to reduce flight test AFCS development 
time and to achieve superior handling qualities characteristics 
through improved design. 

2. Use derivative identification during initial flight tests of 
prototype aircraft to assess and locate any flight dynamic problems. 

3. Use derivative identification to expand the flight envelope of 
present and new aircraft. 

h. Use derivative identification as support in upgrading analytic 
prediction methods so that future, aircraft can be designed more 
accurately. 

To progress toward the above goals, it is essential to pursue applica- 
tions and continue developments and refinements of derivative identification. 
New developments in low speed airspeed instrumentation, automatic instru- 
mentation trimming, and redundant instrumentation are some refinements that 
could be made. The need to separate rotor degrees-of-freedom effects from 
body degrees of freedom requires 9 or 12 DOF identification. The same 
analytic methods developed in this report apply to this problem with minor 
modifications. Continued research in this area can provide for total flight 
dynamic derivative correlation. 

The following specific areas are recommended for further research: 

1. Development and application of a computationally efficient algorithm 
for identification of 9 and 12 DOF models. For the articulated 
rotor helicopter, this allows body, tip-path-plane, and lag DOF. 

The methods developed in this report are applicable to this problem 
with slight modifications. 

2. Incorporation of the fixed interval smoothing solution to yield 
gust and initial condition identification. Improved derivative 
estimates can then be obtained on the second pass, using the 
estimated gusts as additional inputs and using the estimated 
initial conditions. 
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APPENDIX A 


THE BAYESIAN MAXIMUM LIKELIHOOD ESTIMATOR 

The equations used for the MLE are presented in equation (T) ana equation 
(8) of the text. The symbols and terms used in these equations will be 
presented here. 

The system equations for a single maneuver are defined in equation (l) 
and equation (2), the state and parameter vector defined in equation (3) 
and the state veriable form is given by equation (4). For the multi-maneuver 
case the vector function f(X,, X ) is defined in equation (la) for 4 maneuvers. 


f(X, Xp) = 
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f 2 ( £>> V 
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0 
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(la) 


For the 4 maneuver case the definitions given below are made, where the 
subscript pertains to the particular maneuver number. 
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Expansion of the partial derivatives in equation (Ha) results in the 
matricies in the following pages. 
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rH 

rH 

rH 

rH 

rH 

rH 




1 — 1 

CM 

OO 

-4- 

LTN 

VO 




cd 

<J3 


d 


ai 

o 

o 

o 


> -- J 


II 

C 


kl 



a -, j ~ (~rsin<}> + qcos<}>) . tane 

2 

= ^i*cos<j> + qsiuy} . sec 9 
a^g = tanO sin<f> 

a^ = cos<p tanO 

a^ 1 = -qsin<}> -rcos<j> 

a^g = cos<(> 

a 29 = -sin<j» 

a^ = (-rsin<J> + qeos<j> ) /cos0 

a^ 2 = (rcos<|> + qsin<|> ) . tan0 

a^g = sin<j>/cos0 

a^p = cos<f>/cos0 

a^ = -gcos0 

a^ 1 = +gcos0cos<)> 

a^ 2 = -gsin<J> .sin0 

a^^ = -gcos0 . sinij> 

agg = -gcosij) .sin0 
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. seed 

I 

I 


df. 

MATRIX F = 3X , (i represents maneuver number) 
“P 


All elements of F are zero except as defined below. 


F, = 6 V 

41 X 

F, = 6 V 

47 y 

F 4,13 ’ 5V z 

F = 6 V 
52 X 

F 0 = 6 V 
58 y 

F 5.1k = 

Fs- = 6 V 
63 x 

F 69 ■ Sv y 

F 6,15 = w . 

F„« = 6 V 

74 X 

F 7,io ■ 5v y 

F 7,16 * 4V S 

F ft _ = 6 V 
o5 x 

F 8,11 = 5V y 

F 8,17 = 

= 6 V 
96 X 

*9,12 = 4V y 

P 9,18 ' 


^,19 
F 5,20 
F 6, 21 
*7,22 
F 8,23 
* 9,24 


6p 

<5P 

<5p 

6p 

6p 

6p 


F 4,25 " 64 

F 4.31 ■ 5r 

F 5,26 ' 44 

F 5,32 ' 4r 

F 6,27 = 44 

F 6,33 “ 4r 

F 7,28 - 64 

F 7,34 = <r 

F 8,29 * 64 

F 8,35 = 6r 

*9,30 " 

F 9,36 ■ 5r 


F 4,37 = 6B ls 

F 4,4 3 = 6A ls 

F 5,38 = 4B ls 

F 5 ,44 = 6A ls 

F 6,39 = 4B ls 

F 6,i*5 = 6A ls 

F 7,4o * 4B ls 

F T,46 = 6A is 

F 8,4l * 4B ls 

F 8,l+7 = 6A ls 

F 9,42 ' 5B ls 

F 9,48 = 6A ls 


4.49 

? 

5.50 

p. 

'6,51 

p 

'7,52 

'8,53 

? 9,5k 


60 


T 


= 60 


T 


= 60 


T 


= 60. 


= 60. 


= 60. 


F, 


P 4,55 

? 5,56 
'6 , 57 
? 7,58 
8,59 
? 9,6o 


= 60 


= 60 


= 60 


= 60 


= 60 


= 60 


The preceding equations along with equation (l), (2), (3), (7) and (8) of 
the text comprise the Maximum Likelihood Estimator. 
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APPENDIX B 


THE INITIAL DERIVATIVE GUESS USING LEAST SQUARE ESTIMATION 

The MLE as presented in the report requires an initial derivative guess 
and variance in the guess. An accurate guess helps keep the linearizations 
required in the algorithm valid and permits the use of less data to obtain 
the desired accuracy. This latter point is important due to the dimensionally 
large problem for 6 degree-of-freedom identification and reduces computer 
time for a solution. It is for these reasons that considerable effort is 
expended in obtaining a close initial derivative guess. 

It is well known that the LSE provides estimates for the derivatives 
which are biased whenever measurement noise is present. Also the variance in 
the estimate tends to be smaller then is actually the case. In addition the 
presence of process noise causes further error in the variance. 

The use of a digital filter followed by a Kalman filter before using the 
least square estimator minimizes the error due to measurement noise and pro- 
vides a good derivative initial estimate. The variance is still too small 
and must be increased to approximate more accurately the uncertainty in the 
derivative estimate. 

This section presents the development of the LSE which is used to 
initialize the MLE. The development allows processing of more than one 
maneuver and permits the use of a priori derivative information. 

The differential equations of motion including the derivatives to be 
identified are given in equations (l) of this report. When using these 
equations with the LSE only the unknown derivatives are required thus all 
known quantities are placed on the left of the equal sign and the derivative 
terms on the right. The equation containing longitudinal acceleration is 
repeated below in this form. 


V - V + gsinQ 
x xo 



«V 

x 



<5V + 

y 


(lb) 


Let all terms on the left of the equal sign he equal to Z f , the pertur- 
bation state variables and control inputs equal the vector Z and all 
derivatives in equation (lb) equal the vector jd. 

Thus equation (lb) can be written as 

Z’ = P T . Z (2b) 

J J 

where ; j = 1, 2 .... n-maneuvers 


Equation (2b) is in the form required for the LSE. The subscript j 
implies this equation holds for each maneuver j . 

The problem is to determine the parameter vector p which minimizes the 
least square performance index of equation (3b). ~ 


n t f 

J *4 ! (Z ’m. - Z 'j , 2 tf j 1 dt + (Po ' P>V <r c ' P> 

J O J 


(3b) 


where ; 


m. 

J 


W s 


represents the terms to the left of the equa] sign-in equation 

(lb) with the measured data substituted for Band V . 

x 

is defined in equation (2b) 

denotes the maneuver number 

is the noise intensity of Z' 

— m . 


W 


represents the true parameters to he identified 

represents the covariance of the initial parameter estimate p^ 


Equation (3b) takes on a minimum value by setting to zero the gradient 
of J with respect to p and solving the resulting simultaneous equations for 

The resulting solution for p is given in equation (4b) and is repeated 
for each row of derivatives. 

p = (w _1 + E f f z wT 1 zT dt) 1 . (4b) 

° j=l o 3 J J 

wT 1 Z' dt) 

J m. 


(W 


-1 


+ E 


n t 

r z. 

j=l o J 


If no a priori information is available W q approaches <» and W~ =0. 
Also replacing integrals by summations results in the discrete equation (5b) 


z. 

“0 


(k) ) (5b) 


P= ( Z 

j=l 


K 

Z 

k=l 


(k) Z T (k) 7 1 


n 


Z 

2=1 


K 

Z 

k=l 


% 


(k) 



where; k represents the data sample number . 

Equation (5b) was used throughout this study to obtain the least square 
derivative estimate. 


The covariance of the derivative estimate^is just the inverted matrix of 

equation (1+b). Since the noise intensity W. is not usually known before 

J 

hand, this quantity can be determined approximately by the fit error between 

the measured data Z' and p^Z as in equation (6b) for the discrete case. 

— m. _ -m 


W! = 
J 


1 

n 




(Z' (k) - 



(k) 


(6b) 


The covariance of the parameter estimate is then given by equation (7b). 


* n K _ — - 

COV ( p - p) = W' ( z z z (k) z; (k) ) x (7b) 

“ ” J j=l k=l J J 


Equation (7b) was used throughout this study to determine the initial 
derivative variance. 
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APPENDIX C 


THE DATA FILTERS 
(Low Pass 5 Digital and Kalman) 

Low Pass 

The Low Pass filter used in this report is simply a first order transfer 
function converted to discrete form for use on a digital computer. Equation 
(lc) represents the equation used. 

w = 2irf 
c c 

y k = (1 - At . w ) y k-1 + (At . w ) . (lc) 

c c 


where ; 

y represents the output of the filter 
k 

Z represents the input data to the filter 
f is the cutoff frequency 

At is the time between samples 

k represents the data point number 

Digital (Graham) 

This zero phase shift digital filter used is discussed in Reference 8 
in complete detail. The basic weight equation used is given in equation (2c). 

h(k.At) = 2 “ ^ (sin(w^kAt) -+ sin(w c kAt) , ) (2c) 

- (w - w )^ (kAt)^ 
t c 


This equation is evaluated for every data point over an interval 
-K<k>K and then normalized by the sum of all the values. The smoothed output 
ZV for each data point, i, then is obtained by equation (3c). 

+k 

Z'i = £ hfc . Z( K + k ) + i 1 = 1, ..., 2K + 1 (3c) 

“JC 

kT 



c 

~ 2-rrf , f 
c* c 

is the 

cutoff frequency 

w t : 

= 2rt, f t 

is the 

termination frequency 




The Kalman Filter 

State variable definitions 


*L 

= <f> 


*13 = Ma \ 

X 2 

= e 


Xu = bias 

X 3 

= <f< 


5 

x 4 

= v 

X 


X, c = bias 
^15 w 6 

X 5 

n 


^6 = Mas l 

X 6 

= V 

z 


^ = bias 2 

*7 

= P 


ha = t>ias 3 

X 8 

= a 


h.9 = hiaS k 

X 9 

= r 


X = bias 
20 5 

^0 

= bias 

W 1 


x 2i = bias 6 

*11 

= bias 

W 2 


X 22 = bias 7 
X 23 = bias 8 

*12 

= bias 

W 3 


x 2l = bias 9 

kQ 

' 





Measurement definitions Distance to c.g. 


Z 1 

— 

<t> 

m 

i ■ 

X 

- v,$,« 

X 

- 

distance 

to 

C.g 

Z 2 

= 

0 

m 

1 ■ 

y 

- V,3,o 

y 

- 

distance 

to 

c.g 

Z 3 

= 

I'm 

m 

l ■ 

z 

- v,e,o 

z 

- 

distance 

to 

c.g 

Z l+ 

= 

V 

m 

X xl 

- X - Accel 

X 

- 

distance 

to 

c.g 

Z 5 

= 

3 

m 

V 

- X - Accel 

y 

- 

distance 

to 

C.g 

Z 6 

= 

a 

m 

1 zl 

- X - Accel 

z 


distance 

to 

c.g 

Z 7 

= 

P m 

1 x2 

- Y - Accel 

X 

- 

distance 

to 

c.g 

Z 8 

= 

\ 

x y2 

- Y - Accel 

y 

- 

distance 

to 

c.g 

Z 9 

= 

r 

m 

1 z2 

- Y - Accel 

z 

- 

distance 

to 

C.g, 




1 x3 

- Z - Accel 

X 

- 

distance 

to 

C.g 




1 y3 

- Z - Accel 

y 

- 

distance 

to 

C.g, 




1 z3 

- Z - Accel 

z 

- 

distance 

to 

C.g, 


Subscripts 

eg - center of gravity 
I - Instrument 
m - measured 

The equations used in the Kalman filter are based on the general equations 
of motion for a rigid body in flight and the equations relating the state 
variables to the measurements . The equations used to obtain the Kalman 
filtered data will be presented in this section. 
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ine general equations of motion for a rigid body in flight are 


<p = p + (rcos<}> + qsin<f> ) tan6 

6 = qcos<f> - rsin<f> 

'P = (rcosifr + qsin<|>)/cos0 

V = Vr - V q - gsin0 + X 
x y z eg 


V = V p - V r + gsin<J> cos0 + Y 
y z * x 6 eg 


(be) 


V = V q-V p + gcos<j) cos0 + Z 
z x y ^ & eg 


P = L 


eg 


q = M 


eg 


r = N 


eg 


The linear accelerations at the center of gravity can he expressed in 
terms of the accelerometer instrumentation by equation (5c). 


a =a T -|(»)xr + a)x 
eg I L 


(a) X rj] 


(5c) 


Equation (5c) thus becomes 


X 

eg 


eg 


' *1 * [ q ^1 ' r Vl ‘ <q2 * p2)1 I l tsI ’ 1 yl ,r P 1 zl] 


*1 - f ^ 


yi 

“z2 


zl_ 

p Ko + P Q 1 X 2 ~ (P 2 + r2 ) + r <1 ! z2 J 


z 

eg 


Z I - [ P V3 


‘ « **3 + p r X x3 + 9 r V 3 - <p2 + 1 2) ^ 3 ] 


(6c) 


The angular accelerations at the center of gravity are directly expressed 
in terms of the angular accelerometers and are given by equation (7c). 


L = L t 
eg I 

M = M_ 
eg t: 

N = N_ 
eg I 


(7c) 
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The accelerometers will in general contain Doth random errors and bias 
errors. A good estimate of the actual acceleration can be obtained by- 
filtering the accelerometer measurement with a digital filter and subtracting 
a constant bias term. Thus the accelerations as read by the accelerometers 
can be written as 




*E 


bias 


bias 


w; 


bias 

v 

bias 


w c 


(8c) 


where; w^,i-l, 2. ..6 represents process noise and has mean value 
the mean value is obtained by filtering the accelerometers with a digital 
filter. 

Substituting equation (8c) into equation (7c) and equation (6c) and then 
substituting the result into equation ( He ) , yields the general system equa- 
tions for the Kalman filter. Using the state variable notation as defined in 
the list of symbols, the system equations become 


Xj_ = Xj + (X 9 cos \ + X Q sin X^) tan X 2 
X 2 = Xg cos X^ - X^ sin X^ 

X 3 = ^ X 9 COS *1 + X 8 Sin x 2 (9c) 

X H = X 5 X 9 " X 6 X 8 - gsin x 2 - ^lo + hk X zl - ^ V 

+ ( X 8 + X 9 ^ 1 xL ~ X 8 *7 1 yl " X 9 X 7 1 zl + w i ” 1 zl W 5 + \l W 6 

X 5 = X 6 *7 “ X 4 X 9 + gSin C ° S X 2 “ *11 + ^5 1 x2 “ ^3 1 z2 

X 8 1 x2 + ^ + X 9 ^ ^ ~ ^ X 8 1 z2 + W 2 + 1 z2 W H “ 1 x2 M 6 
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Vs - Vu Vs 


x 6 = X l> X 8 ' x 5 V * gcos V oos x 2 - X 12 + 


13 


*8 


* V *9 V 3 - X 8 X 9 V3 + <V + X 8 > Vs + *3 - 

-v 3 * V 

Vk + V 


1 0 W 

y3 i. + 1. 


X 9 “ V 5 + w 6 
Vo * 0 


Vu - 0 


The measurement equations are 


<f> m = <(» + bias^ + n^ 


0 = 0 + bias 0 + n^ 

m 2 2 


= ip + bias- + n_ 

m 3 3 


V = (V 2 T + V' 


m 


T . , 2 T * v 2 T ) 1/2 - 

xl yl zl 


bias ^ + n^ 


-1 


g = tan (V T /V ) + bias,. +n r 
m jrl xl 5 5 

a = tan ^ (V T /V T ) + bias^ + n r 
m zl xl 6 6 


P m = P + bias^ + n^ 


q^ = q + biasg + ng 


r = r + bias rt + n. 
m 9 9 


(10c) 
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where 



c.g. 


+ to X r 


(He) 


Substituting equation (lie) into equation (10c) and rewriting in the 
state variable notations as defined in the list of symbols results in the 
measurement equation (12c). 


Z 1 ‘ h + *16 + "l 

Z 2 - X 2 ♦ hr + n 2 

Z 3 - X 3 + hs + n 3 

Z 4 " ( X 1 , + *8 h - X 9 \ )2 

< x 6 + x i\-h V 2 

Z 5 = tan' 1 (X ? + Xj l x - X. 

Zg = tan' 1 (Xg + Xj l y - Xg 

Z 1 “ + X 22 + n 7 
Zg = x 8 + X 23 + n 8 

Z 9 = X 9 + X 2U + n 9 

Equation (9c) and equation (12c) are the system and measurement equations 
respectively and are rewritten below in vector notation. 


X = f(X, t) + Gw(t) 

(13c) 

Z = h(X, t) + n 

(lUc) 


- (x 5 + x 9 h ' V + 
1/2 + ^9 + 


(12c) 


V /(X 1, + X 8 h ' x 9 V> + x 20 * n 5 


1 )/(*>, + X* 1 - X„ 1 ) + X,, + n 


'Vq X — X 

o z 9 y 


21 
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Equation (l3e) .and (l4c) are precisely in the form for use with the 
extended Kalman filter. 


The extended Kalman filter has the solution presented below. 

• «\i_ i • T , 

(15c) 


(16c) 


The initial conditions are given by 

x(t 0 )=2^ (17c) 

P ( t o ) = P o U8c) 

This equation set is solved by forward processing of the data. At the 
end of the data record an estimate of the bias is available. This bias is then 
used as initial estimate and equations (15c) and (l6c) processed a second 
time. The state estimates are considerably improved on the second pass. 


X = f(X, t) + G w + P 


• -p 

p-Hi -p + p 


3 *lX 


3f 

3X 


- I X 


3X 

T 


. R (Z - H (X) ) 
X 


+ gqgT - p II 


R 


- IX 


-1 3h 
3X 


, T 
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APPENDIX D 


The Unstable Error Variance Equation 


After stability derivatives are identified the linear model can be 
simulated and compared against the flight data. The error between the outputs 
of the derivative model and the test data is governed by an error equation 
and error variance equation which has eigenvalues equal to the eigenvalues of 
the identified derivative. The error variance equation is similar to equation 
(8) of the text with R - 1 set to zero. The equation is rewritten below. 


Term 1 


11 


= A P 
11 11 


+ (A. 


Te: 


•T 

’Pi 




= AP 


11 
Term 

J— 


P ll^ 


Term 2 
+ GQG + 


pi ~lp pp 


Term 3 


( pT a T )T 

l pl V 


+ P 


Pi 


£ 


(Id) 

(2d) 


(to) = Initial State Covariance 


Each term in equation (id) and equation (2d) has the effect indicated 
below. 


Term 1 - The matrix represents the identified stability derivatives. 

The stability of equation (id) and equation (2d) depends on the 
eigenvalues of this matrix. 


Term 2 


This term represents the amount of process noise present and acts 
like a step input to equation (id). 


Term 3 


This term is due to errors in the identified derivatives and acts 
like an input driving equation (id). 


Term 


This term determines the stability of equation (2d). 


Term 5 


This term acts like an input driving term to equation (2d). P 
represents the covariance of the identified derivatives. 


Interpretation of the error variance equation 

The error between the outputs of the identified linear model and the 
test data is governed by an error equation. The actual error can be estab- 
lished by observing the simulated response. The degree to which we can 
believe this error (ie the variance in the error) is governed by equation 
(Id) and (2d). The usefulness of the equation can best be understood by 
example . 
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Suppose the derivatives were identified without error . The terms in 
equation (id) and (2d) would all go to zero except Term 1 and Term 2. If 
the matrix A^ had stable eigenvalues then errors due to ini tial conditions 
should diminish with time and should show perfect correlation of linear model 
output and test data. If process noise is present there can be an error 
after transients die out within an error band governed by equation (id). 

If the matrix A_^ is unstable then the error will diverge and it is 
expected to be within an error band governed by equation (id). 

The importance of this equation is that it tells us on the average 
how accurate we should expect the output of the linear model to match the 
test data. For the case with perfectly identified derivatives and unstable 
eigenvalues the expected error grows with time; due to initial condition 
errors and process noise. If there are errors in the identified derivative 
then this error will be even greater as governed by equation (id) and (2d). 
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Table 1. - Calculation of Measurement Noise Statistics for the Kalman Filter 
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This quantity is Inputed by the Ehgineer when not computed from the Data 














Table 2. - Calculation of Process Noise Statistics for the Kalman Filter 
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Table 3. - Calculation of Initial State and Variance For the Kalman Filter 
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Table 3. - Concluded. 
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Table It. - Description of CH-53A Test Vehicle 
Description of Test Vehicle 


Helicopter Model: Sikorsky Aircraft 

U. S. Navy 

Aircraft BuNo. . 

Number of Engines 
Engine Manufacturer 
Engine Type 


Longitudinal C. G. range 
Aircraft Datum. 


S-65 

CH-53A 

153728 

2 

General Electric 

YT64-GE-12 

T64-GE-413 

328" to 352” 

84" Forward of nose 


Main Rotor: 


Tail Rotor: 


Horizontal 

Stabilizer: 


Weight And 
Inertias : 


Number of Blades 

Diameter 

Airfoil 

Chord 

Total Blade Area 
Disc Area 

Solidity Ratio - effective 
Blade Twist 

Shaft Tilt - longitudinal 


6 

72.225 ft. 

NACA 0011 (mod) 
26.0 in. 

374.95 sq. ft. 
4098.13 sq. ft. 
0.1158 

- 6 ° 

5° 


Number of Blades 

Diameter 

Airfoil 

Chord 

Total Blade Area 
Disc Area 

Solidity Ratio - effective 
Blade Twist 

Pitch Range standard 

extended range 


k 

16 ft. 

NACA 0012 
15*^ in. 

33-12 sq. ft. 
201.1 sq. ft. 
0.1980 

- 8 ° 

-2° to +24° 
-2° to +28° 


Area 

Airfoil Section - tip 
- root 

Aspect Ratio 
Taper Ratio 
Incidence Angle 
Dihedral Angle 


40.0 sq. ft. 
NACA 0012 
NACA 0016 
2.5 
0.6 
+ 3-5° 

+ 5-0° 


W 

Ixx 

iyy 

Izz 

Ixz 


35000 Lbs. 

34250 Ft-Lbs-Sec 
166100 Ft-Lb-Sec 
153600 Ft-Lb-Sec 
Assumed zero 
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Table 5 • 


Description of UH-54B Test Vehicle 
Description of Test Vehicle 


Aircraft Bu. No. 


6401-10005 

Helicopter Model: 

Sikorsky Aircraft 

S-64A 


U. S. Army 

CH-54B 

Aircraft Govt . No . 


69-18467 

Number of Engines 


2 

Engine Manufacturer 


Pratt & Whitney 

Engine Type 


JFTD12A-3 

Longitudinal C.G. range 

328" to 346" 

Aircraft Datum 


34" Forward of nose 

Main Rotor: 

Number of Blades 

6 


Diameter 

72.0 ft. 


Airfoil 

NACA 0010.91 


Chord 

23.65 in. 


Total Blade Area 

352.2 sq. ft. 


Disc Area 

4072.0 sq. ft. 


Solidity Ratio - effective 

0.1021 


Blade Twist 

-7.2° 


Shaft Tilt - longitudinal 

30 


- lateral 

3° 

Tail Rotor: 

Number of Blades 

4 


Diameter 

16 ft. 


Airfoil 

NACA 0012 


Chord 

15*4 in. 


Total Blade Area 

33.12 sq. ft. 


Disc Area 

201.1 sq. ft. 


Solidity Ratio - effective 

0.1980 


Blade Twist 

- 8 ° 


Pitch Range 

-2° to +24° 


Horizontal Stabilizer: Area 

Airfoil Section 
Aspect Ratio 
Taper Ratio 
Incidence Angle 
Dihedral Angle 


25.89 sq. ft. 
NACA 0012 

1.1k 

1.0 

0 ° 


Weight And W 

Inertias : Ixx 

iyy 

Izz 

Ixz 


28000 Lbs 
29329 Ft-Lbs-Sec 
150198 Ft-Lbs-Sec 
130999 Ft-Lbs-Sec 
Assumed zero 
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Table 6. 


Methods Used in Determination of The Most Accurate 
Approach For Derivative Extraction. 


Method 

No. 

ru 1 i. 

-ucx ua riiuci 

For 

LSE 

A Priori 
Derivative 
Guess For 
MLE 

Dat a Fi It er 
For 
MLE 

Derivative 

Extraction 

Method 

1 

First Order 

— 

— 

Least Square 

2 

Digital 

- 

- 

Least Square 

3 

Kalman 

- 


Least Square 

It 

- 

Arbitrary^ ^ 
Arbitrary ^ ^ 
Arbitrary^ 

(2) 

Least Square 

(2) 

Least Square 

(2) 

Least Square 

(2) 

Least Square 

(3) 

Least Square 

(k) 

Least Square 

First Order 

Max. 

Likelihood 

5 


Digital 

Max . 

Likelihood 

6 

— 

Kalman 

Max. 

Likelihood 

7 

First Order 

First Order 

Max. 

Likelihood 

8 

Digital 

Digital 

Max. 

Likelihood 

9 

Kalman 

Kalman 

Max. 

Likelihood 

10 

Kalman 

Digital 

Max. 

Likelihood 

11* 

Kalman 

Kalman 

Max . 

Likelihood 

12* 

Kalman 

Kalman 

Max. 

Likelihood 


(1) See Table 7. for the Arbitrary Guess. 

(2) Variances Multiplied by 100. 

(3) Variance obtained by engineering judgement (See Table 18). 

(h) Variance of Method 11 divided by It. 

* Method 11 and Method 12 were not used in selecting the most accurate 
method, but were used on the CH-5^B data. 
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Table 7* - The Arbitrary Derivative Guess and Variance 

Used to Initialize the Max. Likelihood Method 
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Table 9b. - Identified Control Derivatives from CH-53A Flight Data at 100 Knot Trim Using 
The Ten Methods (Presented in Dimensional Form, British System of Units) 
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Table 10. - RMS Error Between Test Data and Identified Derivative Model 

For the Ten Methods (CH-53A, 100 knots). 



Total RMS Error For All State Variables and Accelerations 












Table 12. - RMS Error Between Test Data and Identified Derivative Model 
For the CH-53A at 150 Knots. 
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Derivatives Identified 
Using the Least Square 
’Presented in Normalize' 







Table lUa - Derivatives Identified from CH-53A Data at 150 knots using the Max. Likelihood Method 
* (Method 9). (Presented in Normalized Form, British System of Units) 


I 
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Table lVb - Derivatives Identified from CH-53A Data at 150 knots using the Max. Likelihood Method 
(Method 9 ). (Presented in Dimensional Form, British System of Units). 
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Table 15. - Derivatives Identified from CH-5UB Data at U 5 knots using the Least Square Method. 
(Method 3) (Using Four 6-second Maneuvers, Presented in Normalized Form). 
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Table 1 6 . - Derivatives Identified from CH-54B Data at 45 knots using the Max. Likelihood Method 
(Method 9) (Using Four 6-second Maneuvers, Presented in Normalized Form). 





Table 18 . - Initial Derivative 1 Sigma Uncertainty Obtained from the Least Square Method 
(Method 3) and Modified using Engineering Judgement. (Used to Initialize the 
MLE of Method 9 , CH-5^B Data, One Maneuver Case). 
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Quantity in Parenthesis obtained by Engineering Judgement. 



Table 19. - Derivatives Identified from CH-^B Data at 45 knots using the Max. Likelihood Method 
(Method ll). (initial Derivative Variance obtained using Engineering Judgement, one 
Maneuver Case, Presented in Normalized Form). 



Table 20a - Derivatives Identified from CH-5hB Data at 4 5 knots using the Max. Likelihood Method 
(Method 12). (Optimum Initial Derivative Variance Used, One Maneuver Case, Presented 
in Normalized Form). 
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Derivatives Identified from CH-54-B Data at 45 knots using the Max. Likelihood Method 
(Method 12). (Optimum Initial Derivative Variance Used, One Maneuver Case, Presented 
in Dimensional Form). 



RIGID BODY UPPER FREQUENCY 



Figure 1. Power Spectrum of Measured Data Showing The Physical Interpretation 
of Measurement Noise Covariance R. 
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Figure 2. Flow Diagram of the Procedure Used to Determine the Most 
Accurate Helicopter Derivative Extraction Method. 












Initializations for the Max. Likelihood Estimator 
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Figure 2. - Concluded 
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Figure 3. CH-53A Test Aircraft 
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■ oti 



a Filtered at 10 HZ 



gure s . - Flight Test Data FROM CH-53A Filtered With A First 

Order Low Pass Filter at 10 HZ. (100 knots, Maneuver 2). 





Filtered 



Continued 






it Test Data Filtered at 10 HZ 



Figure 6 . _ Flight Test Data From CH-53A Filtered With A First 

Order Low Pass Filter At 10 HZ. (100 knots, Maneuver 3). 
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Figure 








Teat Data Filtered 







Flight Test Data From CH-53A Filtered With A First 

Order Low Pass Filter At 10 HZ (100 knots. Maneuver 6). (Data not used 
Identification. ) 









Figure 9 . - Continued 





— Kalman Filtered Flight Data 
Digital Filtered Flight Data At 
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Figure 10 . - Comparison of Flight Test Data From CH-53A using the 

Kalman and Digital Filtering Methods (100 knots. Maneuver l) . 
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Kalman Filtered Flight Data 
Digital Filtered Flight Data At 3 HZ. 
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Figure 13 . - Comparison of Flight Test Data From CH-53A using the 

Kalman and Digital Filtering Methods (100 knots. Maneuver h) . 









Kalman Filtered Flight Data 
Digital Filtered Flight Data At 3 HZ. 
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Figure 13 . - Concluded 
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Figure lit . - Concluded 



.euver 6) . (Data not Used in the 









Figure l6 . - Continued 





Flight Test Data Filtered At 10 HZ 
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Figure 17 . - Flight 


Flight Test Data Filtered At 10 HZ 
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Flight Test Data Filtered At 10 HZ 
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Figure 19 . 




Flight Test Data Filtered At 10 HZ 






Kalman Filtered Flight Data 
Digital Filtered Flight Data At 3 HZ. 
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Figure 20 • - Comparison of Flight Test Data From CH-53A Using The Kalman 

And Digital Filtering Methods (150 knots. Maneuver l) . 





Kalman Filtered Flight Data 
Digital Filtered Flight Data At 3 HZ. 
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Figure 21. - Comparison of Flight Test Data From CH-53A Using The Ka lm a n 

And Digital Filtering Methods (150 knots. Maneuver 2). 
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Figure 22. - Comparison of Flight Test Data From CH-53A Using The Kalman 

And Digital Filtering Methods (150 knots. Maneuver 3). 







Kalman Filtered Flight Data 
Digital Filtered Flight Data At 3 HZ. 
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Figure 23. - Comparison of Flight Test Data From CH-53A Using The Kalman 

And Digital Filtering Methods (150 knots. Maneuver 1*) . 
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Figure 2h. CH-54B Test Aircraft 





Flight Test Data Filtered at 8 HZ 



Figure 25. - Flight Test Data From CH-51+B Filtered With a First Order 

Low Pass Filter at 8 HZ. (1+5 knots, 1 6 sec. Maneuver). 


Flight Test Data Filtered at 8 HZ 
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Figure 25 . - Continued 



Figure 25 - Continued. 
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Figure 25. - Continued 


Flight Test Data Filtered at 8 HZ 



Figure 25. - Concluded 





Kalman Filtered Flight Data 
Digital Filtered Flight Data at 3 HZ 



Continued 






Kalman Filtered Flight Data 
Digital Filtered Flight Data at 3 HZ 
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Kalman Filtered Flight Data 
Digital Filtered Flight Data at 3 HZ. 
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Figure 26. « Continued 


Kalman Filtered Flight Data 
Digital Filtered Flight Data at 3 HZ 



Concluded 
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Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 10) 
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Square Derivative Model (Method 3) 

dkelihood Derivative Model (Method 10 



nnjou* 








.t Data 
ve Model (Me' 
ative Model 






Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 10) 





Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 10) 
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Fi gure 


Filtered Flight Dat 
Square Derivative Mo 
.kelihood Derivative 








Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 10) 
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PERTURBATION OF NONLINEAR MODEL 
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Figure 35. Comparison of the Derivative Mp Identified From Flight Data Using The 
MLE vs. Analytically Predicted Value (CH-53A, 100 Knots). 
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Figure 36 . Roll Damping Derivative Convergence and Uncertainty From the M 
For The Two Different Initial Derivative Starting Values (From 
CH-53A Flight Data, 100 Knots). 
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Figure 37 • Pitch Control Derivative Convergence and Uncertainty From The 
MLE for Two Different Initial Derivative Starting Values (From 
CH-53A Flight Data, 100 Knots). 
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Figure 38. Pitch Damping Derivative Convergence and Uncertainty From The MLE For 
Different Initial Derivative Starting Values (From CH-53A Flight Data, 
100 Knots). 




□ - USING LEAST SQUARE 
IDENTIFIED DERIVATIVES 



Figure 39- Characteristic Roots Obtained From The Identified Derivatives of 
The LSE (Method 3) and The MLE (Method 10 ) (From CH- 5 3A Flight 
Data, 100 Knots, AFCS Off) 
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■ FLIGHT TEST DATA (6 D.O. F. IDENTIFIED) 
<> THEORY (9 DO. F. IDENTIFIED) 
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Figure ho. Comparison of Characteristic Roots Identified From Flight Data 
Analytic Computer Model Roots. (CH-53A, 100 Knots, AFCS OFF). 
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Figure k2. Characteristic Roots Obtained From The Identified Derivatives of 
The LSE (Method 3) and The MLE (Method 9) (From CH-53A Flight 
Data, 150 Knots, AFCS Off). 
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Figure 43. Comparison of Characteristic Roots Identified From Flight Data Vs. Analytic Computer 
Model Roots and Showing Unstable Dutch Roll Roots Obtained by Visually Inspecting 
Independent Test Data. (CH-53A, 150 Knots, AFCS Off). 




• LATERAL DIRECTIONAL MODE ONLY 
0 LONGITUDINAL MODE ONLY 



Figure kk. Characteristic Roots Obtained By Complete Decoupling of Longitudinal 
From the Lateral Identified Derivatives Showing That The Unstable 
Dutch Roll Roots are Caused By Longitudinal-to-Lateral Coupling. 
(From CH-53A Flight Data, 150 Knots, AFCS Off). 
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▲ SIX DEGREE OF FREEDOM COUPLED ROOTS 
A ROOTS WITH N q = 0 



Figure 45. Characteristic Roots Obtained From The Full Identified Derivative 
Array Vs. Roots Obtained When the Derivative Nq is Set To Zero. 
Demonstrating the Destabilizing Effect of Nq. (From CH-53A 
Flight Data, 150 Knots, AFCS Off). 
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Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 9) 




CH-53A Flight Data (150 knots 
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Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 9) 
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Figure ^7 . _ Time History Comparison of Identified Derivative Models Against 

CH-53A Flight Data (150 knots, Maneuver 2). 


Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 9) 




Figure 47. - Continued 
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Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 9) 
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CH-53A Flight Data (150 knots 
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CH-53A Flight Data (150 knots 







Filtered Flight Data 



Time History Comparison of Identified Derivative Models 
Against CH-5^B Flight Data (15 knots, 1 6 sec. Maneuver) 



Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Maxi Likelihood Derivative Model (Method 12) 



Figure 5^. - Continued 
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x. Likelihood Derivative Model (Method 12) 



Figure 50, _ Continued 



Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 12) 
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Figure 50 . _ Continued 


Kalman Filtered Flight Data 

Least Square Derivative Model (Method 3) 

Max. Likelihood Derivative Model (Method 12) 





Figure 50. _ Continued 
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Figure 50 # _ Concluded 
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Figure 51. - Stability Derivative Convergence and 1 Sigma Uncertainty 

for CH-5^B using Method 12. (l6 sec. Maneuver) 
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Figure 51. - Continued 










Figure 51 . - Continued 





Figure 51. - Continued 





34 .000 



Figure 51. - Continued 





Figure 51. - Continued 
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Figure 51 . _ Continued 





Figure 51 . - Continued 





Figure 51 . - Continued 
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Figure 51 . - Continued 
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Figure 51 . - Continued 











SindNI X3I1S 101 Id 


238 


Figure 53. Pilot Input Sequence Recommended. For Test Programs To Generate 
Data For Stability Derivative Identification. 



DERIVATIVE IDENTIFICATION PROCEDURE 
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Figure 5^. Stability Derivative Identification Procedure Found to Be Most 

Accurate For Extracting Derivatives from Helicopter Flight Data. 






